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PREFACE 


This  report  is  one  of  a  series  describing  the  results  of  numerical 
experiments  with  the  two-level  Mintz— Arakawa  model,  which  are  being 
carried  out  as  part  of  the  Rand/ARPA  Climate  Dynamics  Program.  The 
present  work  presents  a  systematic  comparison  of  the  observed  January 
climate  with  a  reference  or  control  integration.  The  identification 
of  the  errors  in  this  simulation  is  a  necessary  preliminary  to  the 
analysis  of  subsequent  experiments  designed  to  explore  the  effect  of 
altered  boundary  conditions.  Other  Rand  publications  related  to  the 
present  report  are  R-877-ARPA,  which  presents  a  detailed  documentation 
of  the  atmospheric  model,  and  R-958-ARPA,  which  compares  the  model  s 
simulated  mean  forcing  fields  with  observation. 
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SUMMARY 


The  mean  global  distributions  of  pressure,  temperature,  wind, 
moisture,  cloudiness,  precipitation,  evaporation,  and  surface  heat 
balance  simulated  for  January  by  the  two-level  Mintz-Arakawa  model 
are  compared  with  the  corresponding  observed  fields.  In  general  the 
simulated  climate  may  be  considered  reasonably  satisfactory,  particu¬ 
larly  in  view  of  the  relative  simplicity  of  the  model.  The  prominent 
errors  are  a  simulated  cloudiness  of  about  half  that  observed  (in  the 
northern  hemisphere) ,  and  a  precipitation  rate  about  twice  that  ob¬ 
served  and  extending  over  too  broad  a  zone  in  the  tropics .  The  cloud¬ 
iness  error  is  evidently  due  to  the  model's  production  of  clouds  only 
during  precipitation,  and  its  failure  to  simulate  nonprecipitating 
cloudiness  at  all.  The  precipitation  error  is  due  to  an  apparent 
simulation  of  excessive  convective  rainfall  and  has  noticeably  affected 
the  heat  and  moisture  balances. 

The  simulated  January  surface  air  temperature  is  higher  than  the 
specified  sea-surface  temperature  over  the  tropical  oceans,  and  has 
resulted  in  a  net  downward  sensible  heat  flux  at  the  surface  between 
about  20°N  and  20° S  in  contrast  with  observation.  The  evaporation, 
which  occurs  almost  exclusively  over  the  oceans,  is  simulated  to  be 
about  50  percent  too  great  in  the  tropics.  These  errors  evidently 
compensate  each  other  in  the  surface  heat  balance  everywhere  except 
at  high  southern  latitudes,  where  the  low  amounts  of  simulated  cloud¬ 
iness  permit  excessive  surface  insolation. 

At  800  mb  the  simulated  relative  humidity  is  approximately  half 
that  observed  between  about  30°N  and  20° S,  while  at  higher  latitudes 
it  exceeds  observation  by  an  average  of  about  15  percent.  At  400  mb 
the  simulated  zonally  averaged  January  temperature  is  approximately 
5  deg  C  above  that  observed  in  the  equatorial  and  tropical  regions; 
at  800  mb  the  simulated  temperature  more  closely  resembles  observation. 
The  meridional  gradients  of  geopotential  height  in  midlatitudes  at 
both  400  mb  and  800  mb  are  somewhat  steeper  than  those  observed,  as  is 
the  meridional  gradient  of  the  400-mb  zonal  mean  temperature.  The 
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associated  maximum  zonal  winds  at  400  mb  and  800  mb  are  about  60  percent 
stronger  than  the  observed  winds,  at  least  in  the  northern  hemisphere. 

These  errors  may  in  turn  be  related  to  the  model’s  tendency  to  sim¬ 
ulate  too  great  a  strength  for  the  quasi-s tat ionary  oceanic  cyclones  of 
middle  and  higher  latitudes,  while  the  intensity  of  the  individual 
transient  waves  is  underestimated.  The  associated  midlatitude  Ferrel 
cell  in  the  mean  meridional  circulation  is  therefore  both  too  weak  and 
too  narrow.  The  subtropical  oceanic  anticyclones  are  more  realisti¬ 
cally  simulated,  and  at  least  in  the  northern  (winter)  hemisphere,  the 
strength  of  the  associated  Hadley  circulation  resembles  that  derived 
from  observations. 

A  more  extended  summary  and  a  discussion  of  the  research  under  way 
on  the  major  sources  of  error  in  the  model  are  given  in  Sec.  VII. 
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I .  INTRODUCTION 


Numerical  models  of  the  atmospheric  circulation  have  been  used  to 
study  the  dynamics  of  climate  since  the  pioneering  calculation  of 
Phillips  (1956) .  The  groups  with  the  most  extensive  experience  in  this 
work  are  at  the  Geophysical  Fluid  Dynamics  Laboratory  of  NOAA,  at  the 
National  Center  for  Atmospheric  Research  (NCAR) ,  and  in  the  Department 
of  Meteorology  at  UCLA.  Starting  from  the  relatively  simple  hemispheric 
model  integrations  reported  by  Smagorinsky  (1963) ,  the  Princeton  group 
has  developed  a  hierarchy  of  models  of  the  general  circulation  with 
increasing  physical  and  numerical  resolution.  In  the  most  recent  (nine- 
level)  formulation  published,  Holloway  and  Manabe  (1971)  have  discussed 
the  comparison  of  the  global  simulation  of  selected  climatic  elements 
with  observation.  The  results  of  the  application  of  the  most  recent 
(six-level)  version  of  the  NCAR  model  to  the  simulation  of  the  global 
climate  have  been  reported  by  Kasahara  and  Washington  (1971) .  These 
comparisons  reveal  a  fair  degree  of  correspondence  between  the  simulated 
and  observed  climate,  and  are  an  important  measure  of  the  models* 
fidelity. 

The  two-level  general-circulation  model  developed  at  UCLA  by  Mintz 
and  Arakawa  beginning  in  1961  has  also  received  extensive  use.  Al¬ 
though  this  model  has  fewer  levels  than  the  NOAA  and  NCAR  models,  it 
contains  a  comparable  degree  of  physical  sophistication,  particularly 
in  the  treatment  of  radiation,  cloudiness,  and  convection  (Arakawa, 
Katayama,  and  Mintz,  1969).  A  preliminary  view  of  the  simulations  per¬ 
formed  with  an  early  version  of  this  model  has  been  given  by  Mintz 
(1965,  1968). 

It  is  the  purpose  of  the  present  report  to  present  a  comprehensive 
comparison  with  observation  of  the  global  climatic  simulations  made 
with  a  recent  formulation  of  the  two-level  Mintz-Arakawa  model  (Gates 
et  al.,  1971).  The  performance  of  this  model  may  then  be  evaluated 
alongside  those  of  the  other  models  cited.  Because  of  its  economy  and 
availability,  this  model  has  been  selected  for  use  in  a  program  of  ex¬ 
perimental  integrations  whose  aim  is  the  study  of  the  response  of  global 
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climate  to  selected  perturbations  and  to  changes  in  the  model’s  physical 
boundary  conditions.  To  evaluate  the  significance  of  such  simulations 
of  climate  modification  properly,  the  model’s  errors  in  specifying  a 
reference  or  control  climate  must  first  be  known. 
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II ,  SUMMARY  OF  THE  MODEL 

The  Mintz-Arakawa  two- level,  atmospheric  general-circulation  model 
employs  the  primitive  (hydrostatic)  equations  at  two  tropospheric  levels 
surrounding  the  globe.  This  model  was  summarized  in  its  essentials  by 
Arakawa,  Katayama,  and  Mintz  (1969)  ,  and  a  comprehensive  description 
and  documentation  of  both  the  model  and  its  numerical  solution  were  re¬ 
cently  prepared  by  Gates  et  al.  (1971).  In  broad  outline,  it  predicts 
the  horizontal  wind  velocity,  temperature,  and  pressure  (or  geopotential) 
at  each  of  two  levels,  which  approximate  the  400-mb  and  80 Q-mb  surfaces, 
with  the  water-vapor  mixing  ratio  predicted  only  at  the  lower  level.  The 
vertical  motion  as  given  by  a  (see  next  page)  is  determined  at  a  central 
level  (approximately  600  mb) ,  and  is  defined  to  be  zero  at  the  top  of 
the  model  atmosphere  (200  mb)  as  well  as  at  the  earth’s  surface  (where 
the  geopotential  is  assumed  known) .  Identification  of  the  model’s  levels 
is  given  in  Fig.  1. 
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Fig.  1  —  Schematic  representation  of  the  model’s  vertical  structure. 
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The  diabatic  heating  in  the  model  is  composed  of  the  net  absorbed 
short-  and  long-wave  radiation,  the  release  of  latent  heat  during  con¬ 
densation,  and  the  sensible  heating  from  the  surface.  The  radiative 
heating  is  itself  affected  by  the  predicted  distribution  of  water  vapor 
and  clouds >  with  the  latter  determined  as  a  result  of  both  large-scale 
and  convective  processes. 

The  moisture  source  is  surface  evaporation  and  the  sink  is  the 

precipitation  released  in  both  large-scale  and  convective  condensation. 

There  is  a  momentum  transfer  by  skin  friction  at  the  earth’s  surface, 

and  an  internal  frictional  momentum  transfer  between  the  two  levels, 

which  depends  on  the  vertical  wind  shear;  there  is  no  lateral  eddy 

♦ 

viscosity  or  diffusion.  In  addition  to  the  requirement  0  -  0  at  the 
earth's  surface,  a  surface  thermal  boundary  condition  of  zero  net  heat 
flux  is  employed  over  land  (and  ice-covered)  surfaces,  and  is  used  as 
the  basis  for  the  determination  of  the  ground  temperature  itself.  Over 
the  oceans  the  sea-surface  temperature  is  prescribed  throughout  the 
simulation.  In  the  present  January  simulation,  the  topography  (land 
elevation,  ocean  surface  temperature,  and  the  locations  of  land  and 
sea  ice)  is  as  given  by  Gates  et  al.  (1971)  in  the  documentation  of 
the  model. 

It  is  here  sufficient  to  state  the  model’s  governing  equations 
for  the  primary  dependent  variables,  and  to  summarize  the  determination 
of*  those  elements  of  the  momentum,  hydrologic,  and  heat  balances  whose 
simulated  distributions  are  shown  later.  We  first  define  the  a-coordi- 
nate  system  by 


o  -  (p  -  pt)/tt  (1) 

where  p  is  pressure,  p^  an  assumed  tropopause  pressure  (200  mb) ,  and 

tt  *  p  -  p-  is  a  measure  of  the  variable  surface  pressure,  p  .  In 
s  1  ® 

this  system,  the  earth’s  surface  (p  =  p  )  is  everywhere  given  by  0  =  1, 

s 

and  the  top  of  the  model  atmosphere  (p  =  p^,)  is  given  by  O  =  0 .  The 
model’s  equations  of  motion  and  thermodynamic  energy  may  be  written 

(ttv)  +  V(TTV,V)  ±  ir a2V2  +  irfkxV  +  trVcj)  +  cnraVir  -  irF 


(2) 
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|^-  (ttT)  +  V  •  (ttTV)  ±  (p/po)lcira202  -  (irac/cp)  (|£  +  V  *  Vir)  =  TtH/cp  (3) 

The  notation  V(tTV,V)  is  to  be  interpreted  as  V(V*7TV)  +  (t TV*V)V.  Here 
V  is  the  horizontal  velocity,  <j>  the  geopotential,  a  the  specific  volume, 
f  the  Coriolis  parameter,  T  the  temperature,  and  0  the  potential  tem¬ 
perature,  and  F  and  H  are  the  frictional  force  and  diabatic  heating  rate 
per  unit  mass.  When  the  equations  are  applied  at  the  model’s  upper 
level  (i.e.,  level  1,  or  O  «  <7^  »  1/4),  the  upper  signs  (plus)  of  the 
third  terms  of  Eqs .  (2)  and  (3)  are  used;  when  they  are  applied  to  the 
lower  level  (level  3,  or  a  *  =  3/4),  the  lower  signs  (minus)  are  used. 

At  these  levels  the  variables  V,  <f>,  a,  F,  T,  p,  and  H  assume  the  sub - 
scripts  1  and  3,  respectively.  In  (2)  and  (3)  the  velocity  and  the 
potential  temperature  02  at  the  midlevel  (i.e.,  level  2,  or  a  »  -  1/2) 

are  taken  to  be  the  average  of  the  values  at  levels  1  and  3.  The  geo¬ 
potential,  (j>,  is  found  by  integration  of  the  hydrostatic  equation, 

9<f)/9cf  =  -ira,  with  the  assumption  that  the  potential  temperature, 

0  *  T(pg/p)K,  is  linear  in  p^-space,  where  K  =  0.286  and  p^  =  1000  mb. 

At  the  lower  level,  the  water-vapor  continuity  equation  is  applied 
in  the  form 

It  (irq3)  +  V  •  irq3V3>5  =  2g(E  -  C)  (4) 

where  q^  is  the  mixing  ratio  at  level  3,  ^  is  the  wind  velocity  ex¬ 

trapolated  to  the  (fictitious)  surface  Cf  =*  7/8  (level  3.5),  and  E  -  C 
is  the  difference  between  the  surface  evaporation  rate,  E,  and  the 
condensation  rate,  C.  The  precipitation  rate  is  assumed  to  be  equal  to 
C,  with  no  allowance  made  for  either  the  liquid  water  content  of  clouds 
or  for  evaporation  from  falling  rain  or  snow. 

At  the  earth’s  surface  the  rate  of  change  of  the  surface  pressure 
as  obtained  from  the  continuity  equation  may  be  written  in  the  form 

=  -V  •  1^2 


9tt 

9t 


(5) 
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while  at  level  2  the  pressure  tendency  may  be  written  in  terms  of  0 
in  the  form 

a 2  =  -(4TT)-1  V  •  TT^  -  V3)  (6) 

The  finite-difference  approximations  used  to  solve  the  system  (2) 
to  (6)  are  due  to  Arakawa,  and  are  described  in  detail  elsewhere  (Gates 
et  al. ,  1971;  see  also  Arakawa,  1972).  After  the  introduction  of  ap¬ 
propriate  scale  factors  for  mapping  the  sphere,  the  general  computa¬ 
tional  philosophy  is  to  write  the  difference  equations  in  such  a  way 
that  the  model’s  total  mass  and  vorticity  are  exactly  conserved,  the 
total  energy  and  momentum  are  conserved  in  the  absence  of  sources  and 
sinks,  and  the  mean  potential  temperature  (and  its  square)  are  very 
nearly  conserved.  The  computational  grids  consist  of  two  interlocking 
set^  of  points  with  a  resolution  of  4  deg  latitude  and  5  deg  longitude. 
One  grid,  on  which  7T,  T,  and  q  are  carried,  includes  the  points  at  ±2, 

±6,  ±10,  ...  deg  latitude  and  0,  ±5,  ±10,  ...  deg  longitude;  the  other 
grid,  on  which  the  wind  velocity  components,  u  and  v,  are  carried,  in¬ 
cludes  the  points  at  0,  ±4,  ±8,  ...  deg  latitude  and  ±2.5,  ±7.5,  ±12.5, 
...  deg  longitude.  The  finite-difference  analogs  are  formed  by  appro¬ 
priate  averages  and  interpolations  in  the  various  centered  and  uncentered 
space-difference  schemes.  The  time  integration  proceeds  in  steps  of  6 
minutes  by  means  of  a  version  of  the  backward  differencing  scheme  due 
to  Matsuno,  which  introduces  a  damping  of  the  solutions  of  high  fre¬ 
quency  in  time. 

The  source  terms,  F,  H,  and  (E  -  C)  in  (2)  to  (4)  are  calculated 
every  fifth  time  step,  at  which  time  an  areal  smoothing  of  the  local 
lapse  rate  is  performed.  The  frictional  force,  F,  consists  of  an  in¬ 
ternal  stress  at  levels  1  and  3  (proportional  to  the  vertical  wind 
shear),  and  a  surface  stress  or  drag  (dependent  upon  the  surface  wind, 
taken  as  0.7  of  the  wind  extrapolated  to  level  4).  The  surface  drag 
acts  on  the  wind  at  level  3.  The  heating  rate,  H,  at  level  1  contains 
the  absorption  of  solar  radiation  in  the  upper  layer  (dependent  upon 
humidity  and  clouds),  the  net  long-wave  radiation  lost  from  the  layer 
(dependent  upon  temperature,  humidity,  and  clouds),  and  the  warming  of 
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the  layer  by  the  latent  heat  released  in  the  parameterized  midlevel 
and  penetrating  convection.  The  heating  rate  at  level  3  contains  the 
corresponding  quantities  for  the  lower  layer,  and  in  addition  includes 
both  the  surface  sensible-heat  flux  (dependent  upon  the  ground-to-air 
temperature  difference  and  the  surface  wind)  and  the  heating  due  to  the 
latent  heat  released  in  large-scale  condensation. 

In  the  model  the  convective  precipitation  depends  upon  a  parameter¬ 
ized  convective  adjustment  process  governed  by  the  vertical  stability, 
while  the  large-scale  precipitation  depends  upon  the  apparent  super¬ 
saturation  (primarily  a  result  of  large-scale  vertical  motion) ,  Mid¬ 
level  clouds  in  the  model  always  accompany  large-scale  precipitation, 
but  the  convective  clouds  are  separately  parameterized  in  terms  of  the 
relative  humidity  at  level  3.  An  analysis  of  the  mean  January  forcing 
fields  given  by  the  net  surface  frictional  force  (surface  stress) ,  net 
diabatic  heating  (H^  +  H^) ,  and  net  moisture  addition  rate  (E  -  C)  is 
given  elsewhere  (Gates,  1972). 
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III.  THE  CONTROL  INTEGRATION 


In  order  to  provide  a  reference  case  against  which  to  judge  the 
significance  of  later  integrations  in  which  deliberate  experimental 
changes  are  made  in  the  model*  in  its  associated  boundary  conditions, 
or  in  both,  a  30-day  period  has  here  been  adopted  as  a  control  simula¬ 
tion  for  comparison  with  January  observations.  In  terms  of  the  calendar 
position  of  the  sun,  this  period  extends  from  29  December  to  27  January; 
in  model  time  this  period  is  from  day  370  to  day  400  of  a  longer  simula¬ 
tion  started  from  an  initial  state  of  rest  and  isothermality  by  Mintz 
and  Arakawa. 

Of  the  many  constants  and  parameters  specified  in  the  model,  atten¬ 
tion  is  here  called  only  to  the  prescribed  boundary  conditions.  As  was 
previously  noted,  the  top  of  the  model  atmosphere  is  taken  as  the  200-mb 
isobaric  surface.  Here  a  =  0  is  required,  and  the  incoming  solar  radia¬ 
tion  corresponding  to  a  solar  constant  of  2  ly  min  ^  is  partitioned  into 
65  percent  subject  to  atmospheric  scattering  and  35  percent  subject  to 
absorption.  At  the  earth’s  surface  the  distributions  of  land  elevation 
and  sea-surface  temperature  are  assumed  to  be  given  as  shown  in  Figs .  2 
and  3,  which  also  serve  to  document  the  continental  outlines  resolved 
by  the  grid.  The  surface  elevations  shown  in  Fig.  2  are  based  on  the 
data  given  by  Berkofsky  and  Bertoni  (1955),  and  the  ocean  surface  tem¬ 
peratures  shown  in  Fig.  3  were  found  from  the  average  annual  data  given 
by  Dietrich  (1963).  In  these  figures  the  grid  points  (of  the  TT-  or 
primary  grid)  are  shown  over  the  oceans ,  and  the  overprinted  symbol  I 
denotes  those  points  at  which  either  ice-covered  land  (Fig.  2)  or  sea 
ice  (Fig.  3)  was  assumed  to  be  present  throughout  the  simulation.  The 
albedo  of  the  surface  is  prescribed  as  0.14  for  bare  land  and  0.07  for 

the  oceans;  the  albedo  of  snow-  or  ice-covered  surfaces  is  designed  to 

* 

vary  between  0.4  and  0.9,  depending  upon  the  latitude  and  season. 

An  inadvertent  program  error  (first  pointed  out  by  E.  S.  Batten) 
has  caused  the  systematic  underestimate  of  the  albedo  of  snow  and  ice 
by  about  0.4.  This  error  is  probably  unimportant  in  the  northern  hem¬ 
isphere  (there  being  little  solar  radiation  north  of  about  70°N  in 
January) ,  but  has  introduced  some  error  over  Antarctica  as  noted  in 
some  of  the  figures,  starting  with  Fig.  6. 
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For  further  details  of  the  model’s  configuration  in  the  control  inte¬ 
gration,  as  well  as  a  complete  and  annotated  listing  of  the  program, 
reference  should  be  made  to  the  model  documentation  (Gates  et  al . , 
1971). 


IV.  THE  SIMULATED  AVERAGE  PRESSURE,  TEMPERATURE,  AND  WIND 


Of  the  many  quantities  computed  in  the  model,  pressure,  temperature, 
and  wind  exhibit  behavior  of  particular  interest;  not  only  are  these  the 
model’s  primary  dependent  variables,  but  they  are  the  meteorological 
variables  most  easily  compared  with  observation.  The  distributions  of 
humidity,  cloudiness,  precipitation,  evaporation,  and  elements  of  the 
surface  heat  balance  are  considered  in  subsequent  sections.  Except  as 
noted,  all  of  the  average  distributions  shown  here  were  determined  from 
data  every  6  hours  during  the  30-day  control  simulation.  The  observed 
distributions  are  from  various  sources  as  compiled  by  Schutz  and  Gates 
(1971) . 

SEA-LEVEL  PRESSURE 

The  global  distribution  of  the  average  January  sea-level  pressure 
simulated  in  the  control  experiment  is  shown  in  Pig.  4,  together  with 
the  observed  average  January  pattern.  In  general  the  model  has  correctly 
depicted  the  principal  centers  of  low  and  high  pressure  over  the  globe. 
The  mean  anticyclones  in  the  subtropics  of  the  southern  (summer)  hemis¬ 
phere  in  particular  are  correctly  simulated  just  to  the  west  of  the 
continents  and  with  approximately  the  correct  maximum  pressure.  The 
Aleutian  low  in  the  northern  hemisphere  is  simulated  about  15  mb  too 
low,  although  its  position  is  reasonably  accurate.  The  Icelandic  low, 
on  the  other  hand,  is  simulated  some  1000  km  southwest  of  its  observed 
average  position,  as  well  as  being  about  15  mb  too  deep. 

The  global  distribution  of  the  difference  between  the  simulated  and 
observed  sea-level  pressure  is  shown  in  Fig.  5.  In  addition  to  the  in¬ 
tensity  errors  in  the  Aleutian  and  Icelandic  lows  noted  above,  large 
differences  in  sea- level  pressure  occur  over  the  continents  of  the 
northern  hemisphere,  where  errors  of  10  to  20  mb  are  common.  It  is 
significant  to  note  that  these  errors  resemble,  in  general  scale  and 
magnitude,  that  found  between  successive  Januaries  in  a  three-year  sim¬ 
ulation  (Mintz  et  al.,  1972),  although  the  simulated  Januaries  all 
show  similar  intensity  errors  in  comparison  with  observation.  This 
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serves  to  emphasize  that  difference  maps  such  as  Fig,  5  strongly  re¬ 
flect  the  errors  of  position  as  well  as  of  magnitude  of  the  quasi¬ 
stationary  (monthly  average)  circulations. 

The  most  serious  error  is  a  systematic  overprediction  of  the  in¬ 
tensity  of  the  semi-permanent  surface  low-pressure  systems,  which 
steepens  the  pressure  gradients  in  Fig.  4a  as  compared  with  those  of 
Fig.  4b  in  the  northern  middle  latitudes.  This  error  is  particularly 
noticeable  in  the  zonally  averaged  sea- level  pressure  given  in  Fig.  6, 
where  the  average  simulated  pressure  at  55 °N  is  about  10  mb  too  low. 

The  only  other  significant  discrepancies  occur  south  of  60° S.  Here 
simulation  exceeds  observation  by  several  millibars  between  about  55°S 
and  70°S,  and  the  sea-level  meridional  pressure  gradient  has  the  wrong 
sign  between  70°S  and  86°S.  This  latter  feature  is  attributed  to  the 
erroneous  surface  heating  over  the  Antarctic  sea  ice  and  Antarctic  con¬ 
tinent  due  to  an  inadvertent  underestimate  of  the  albedo  of  snow  and 
ice.  Over  the  bulk  of  the  atmosphere,  however,  the  simulated  latitudi¬ 
nal  pressure  profile  is  in  good  agreement  with  the  observed.  Since 
the  model  is  constrained  to  have  a  global  average  surface  pressure 
equal  to  that  observed  (984  mb;  see  Gates  et  al.,  1971),  the  difference 
between  the  simulated  global  average  sea-level  pressure  (1013.3  mb)  and 
the  observed  value  (1011.1  mb)  reflects  only  errors  in  reduction  to 
sea  level. 

This  comparison  of  predicted  and  observed  sea-level  pressure  ma^ 
be  considered  to  update  that  given  by  Mintz  (1965,  1968)  for  an  earlier 
version  of  this  model  (which  used  a  7-deg  latitude  by  9-deg  longitude 
grid,  did  not  predict  the  water  vapor  or  cloudiness,  and  had  simpler 
treatments  of  convection  and  radiative  heating,  among  other  differences). 
The  present  model  version  may  be  considered  to  produce  a  more  satis¬ 
factory  overall  simulation.  Its  use  of  a  smaller  grid  size  without 
lateral  diffusion  has  reduced  the  truncation  error  of  the  middle-latitude 
transient  cyclones  in  particular.  Significant  errors  remain,  however, 


This  conjecture  was  confirmed  in  a  subsequent  January  simulation 
in  which  the  surface  albedo  for  snow  and  ice  was  corrected  (Mintz, 
personal  communication);  see  Fig.  6. 


Fig.  4a  -  The  simulated  sea- level  pressure  for  January,  with  isolines  every  5  mb  and  the 
1000— mb  isobar  dashed.  The  symbols  H  and  L  denote  local  maxima  and  minima. 
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Fig.  4b  The  observed  sea-level  pressure  with  isolines  every  5  mb  and  the  1000-mb  isobar  dashed. 
From  Schutz  and  Gates  (1971) ,  based  upon  data  of  Crutcher  and  Meserve  (1970)  and 
Taljaard  et  al.  (1969). 


-level  pressure  (mb) 


6  —  The  zonal  average  sea-level  pressure  as  simulated  (full  line)  and  as  observed  (dashed  line) , 
based  on  the  data  of  Fig.  4.  The  dotted  curve  at  high  southern  latitudes  in  this  and  later 
figures  is  the  simulation  with  a  corrected  albedo  for  snow  and  ice,  as  obtained  from  a  sub¬ 
sequent  January  integration  (Mintz,  private  communication);  see  footnote,  p.  8. 
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in  the  amplitude  (and  location)  of  the  time-averaged  or  quasi-stationary 
systems  (see  below,  and  also  Gates,  1972), 

GEOPOTENTIAL  HEIGHTS  AT  8QQ  MB  AND  400  MB 

The  average  January  geopotential  height  simulated  at  800  mb  and 

400  mb  are  shown  in  Figs.  7  and  10,  along  with  the  observed  global 
& 

pattern.  As  was  the  case  with  the  sea-level  pressure,  the  model  has 
reproduced  the  general  pattern  of  geopotential  height  in  the  tropo¬ 
sphere  with  reasonable  accuracy,  although  the  simulated  mean  troughs 
and  ridges  are  too  large  in  amplitude.  In  the  southern  hemisphere  the 
simulated  geopotential  height  gradient  approximately  agrees  with  obser¬ 
vation  at  both  800  mb  and  400  mb,  as  do  the  pole-to-equator  height 
differences  in  the  northern  hemisphere.  The  800-mb  and  the  400-mb  sim¬ 
ulated  height  gradients  in  the  northern  middle  latitudes,  however,  are 
approximately  twice  as  great  as  those  observed.  Together  with  the  sea- 
level  pressure  (Fig.  4),  these  data  show  the  model’s  ability  to  develop 
the  correct  vertical  tilt  and  thermal  wind  in  the  middle- latitude 
systems . 

The  error  fields  of  the  simulated  geopotential  heights  at  800  mb 
and  400  mb  are  shown  in  Figs.  8  and  11.  We  note  that  the  simulation 
errors  or  discrepancies  are  positive  nearly  everywhere  and  are  several 
hundred  meters  in  magnitude  at  both  levels.  This  error  is  even  more 
clearly  shown  in  the  latitudinal  variations  of  the  zonal  averages  in 
Figs.  9  and  12;  here  the  simulated  800-mb  and  400-mb  heights  are  sys¬ 
tematically  higher  than  the  observed  values  at  all  latitudes.  On  a 
global  average,  the  mean  simulated  heights  at  800  mb  and  400  mb  are 
2268  m  and  7600  m,  which  exceed  the  observed  mean  heights  of  these 
surfaces  by  330  m  and  322  m,  respectively.  The  source  of  this  dis¬ 
crepancy  is  not  the  interpolation  errors  involved  in  determining  the 
simulated  800-mb  and  400-mb  geopotential  heights,  nor  is  it  that  in¬ 
volved  in  finding  the  observed  heights  (see  Schutz  and  Gates,  1971). 
Rather,  it  lies  in  the  determination  of  <J>^  and  themselves  as  devia¬ 
tions  from  the  height  of  the  center  of  mass  of  the  model  atmosphere 


Figures  in  some  sections  have  been  sequenced  to  facilitate  com¬ 
parison  rather  than  in  the  order  cited. 
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(see  Gates  et  al.,  1971)  •  This  height  has  been  selected  in  such  a  way 
that  when  and  are  equal  to  the  observed  temperatures,  the  verti¬ 
cal  integral  of  the  model’s  potential  and  internal  energy  approximates 
the  vertical  integral  of  the  potential  and  internal  energy  of  the  actual 
atmosphere •  The  model’s  center  of  mass  thus  turns  out  to  be  about  326  m 

higher  than  that  of  the  real  atmosphere,  and  accounts  for  the  systematic 

* 

error  seen  in  Figs.  7  to  12. 

SURFACE  AIR  TEMPERATURE 

The  average  January  surface  air  temperature  simulated  in  the  con¬ 
trol  experiment  is  shown  in  Fig.  13,  along  with  the  observed  January 
pattern.  In  broad  terms  the  model  has  correctly  simulated  many  observed 
features  of  the  surface  air  temperature,  of  which  the  most  prominent 
are  the  cold  air  over  the  interior  of  the  continents  of  the  northern 
(winter)  hemisphere,  the  relatively  warm  air  over  the  southern  hemis¬ 
phere  continents,  and  the  location  of  the  freezing  isotherm  (dashed 
line  in  Fig.  13).  Except  in  the  Antarctic  region,  where  there  was  an 
underestimate  of  the  surface  albedo,  the  observed  gradient  of  surface 
air  temperature  is  also  reasonably  well  portrayed  in  the  simulation. 

Over  the  oceans  the  prescription  of  a  fixed  sea-surface  tempera¬ 
ture  has  constrained  the  model’s  simulation  of  the  surface  air  tem¬ 
perature  to  be  within  about  3  deg  C  of  the  prescribed  water  temperature, 
as  may  be  seen  in  the  global  error  of  Fig.  14.  Over  most  of  the  tropical 
oceans,  however,  the  simulated  surface  air  temperature  is  higher  than  the 
sea-surface  temperature  itself  (see  Fig.  3).  This  feature  has  affected 
the  surface  heating  as  is  discussed  in  Sec.  VI  below.  Over  the  con¬ 
tinents,  where  a  surface  heat  balance  was  assumed  in  the  model,  even 
larger  errors  occur.  In  general,  the  land  regions  of  middle  and  high 
latitudes  are  simulated  to  be  too  warm  by  about  10  deg,  while  in  the 
tropics  the  land  areas  are  characteristically  simulated  to  be  several 
degrees  too  cold.  This  is  probably  the  result  of  the  model’s  inability 
to  accurately  portray  either  a  strong  surface  inversion  or  a  steep  low- 
level  lapse  rate.  The  model’s  failure  to  develop  a  strong  surface 

I  am  indebted  to  Yale  Mintz  and  Akio  Arakawa  for  pointing  out 
this  property  in  their  model. 


Fig.  7b  —  The  observed  geopotential  height  at  800  mb  for  January  with  isolines  every  50  m  and 
the  2000-m  contour  dashed.  From  Schutz  and  Gates  (1971),  based  on  data  of  Crutcher 
and  Meserve  (1970)  and  Taljaard  et  al.  (1969). 


Geopotential  height  at  800  mb  ( 


Fig.  9  —  The  zonal  average  800-mb  geopotential  height  as  simulated  (full  line)  and  as  observed 

(dashed  line),  based  on  the  data  of  Fig.  7.  See  Fig.  6  for  dotted— line  identification. 


Fig.  10a  --  The  simulated  geopotential  height  at  400  mb  for  January,  with  isolines  every  100  m 
and  the  7000-m  contour  dashed. 


ie  10b  —  The  observed  geopotential  height  at  400  mb  for  January ,  with  isolines  every  100  m 
8'  and  ?he  7000-m  eontour  dashed.  From  Schutz  and  Gates  (1971).  based  on  data  of 

Crutcher  and  Meserve  (1970)  and  Taljaard  et  al.  (1969). 


i 

ho 

Ln 

I 


Fig.  12  —  The  zonal  average  400-mb  geopotential  height  as  simulated  (full  line)  and  as  observed 
(dashed  line),  based  on  the  data  of  Fig.  10. 


4- 


.  13b  —  The  observed  surface  air  temperature  for  January,  with  isolines  at  5-deg  C  intervals 
and  the  0  deg  C  isotherm  dashed.  From  Schutz  and  Gates  (1971),  based  on  data  of 
Crutcher  and  Meserve  (1970)  and  Taljaard  et  al.  (1969). 
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Fig.  15  —  The  zonal  average  surface  air  temperature  as  simulated  (full  line)  and  as  observed  (dashed 
line),  based  on  the  data  of  Fig.  13.  Also  shown  are  the  zonal  averages  of  the  prescribed 
sea-surface  temperature  (as  used  in  the  simulation)  and  of  the  observed  sea-surface  tempera¬ 
ture  for  January  from  Schutz  and  Gates  (1971),  based  on  data  of  Washington  and  Thiel  (1970). 
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inversion  during  the  polar  night  is  also  the  cause  of  the  high  simulated 

surface  air  temperatures  over  the  Arctic.  The  large  positive  errors,  in 

the  range  of  10  to  20  deg  C,  in  the  surface  air  temperature  in  the  south 

polar  region  are  mainly  due  to  the  model’s  use  of  an  erroneously  low 

* 

albedo  for  snow  and  ice. 

On  the  basis  of  the  zonal  averages  shown  in  Fig.  15,  the  model  is 
seen  to  simulate  a  surface  air  temperature  slightly  too  high  at  most 
latitudes,  leading  to  a  global  average  which  is  1.2  deg  C  above  the 
observed  global  mean  value  of  12.2  deg  C.  Only  between  20°S  and  65°S 
does  the  model  simulate  too  low  a  zonal  average  surface  air  temperature, 
which  may  be  attributed  to  the  prescription  of  too  low  an  ocean-surface 
temperature.  This  may  be  seen  in  Fig.  15,  where  the  (fixed)  sea-surface 
temperature  used  in  the  simulation  is  shown,  along  with  the  zonal  aver¬ 
age  of  the  observed  values  for  January;  where  the  prescribed  ocean 
temperature  is  too  high,  the  simulated  surface  air  temperature  is  also 
too  high.  Use  of  the  observed  sea-surface  temperature  data  for  January 
(rather  than  the  annual  average  used  here)  would  remove  this  error  and 
might  improve  the  simulation  in  the  northern  hemisphere  as  well. 

TEMPERATURE  AT  800  MB  AND  400  MB 

The  simulated  and  observed  average  January  temperatures  for  800  mb 
and  400  mb  are  shown  in  Figs.  16  and  19,  respectively.  The  major  large- 
scale  features  of  the  temperature  field  are  seen  to  be  reasonably  well 
simulated  by  the  model  at  both  levels.  The  principal  axes  of  the  mean 
January  cold  and  warm  air  tongues  are  in  approximately  the  observed 
positions,  with  the  longitude  of  the  maximum  temperature  gradient  cor¬ 
rectly  simulated  off  the  east  coasts  of  North  America  and  Asia.  In 
general,  however,  the  model  has  given  too  large  an  amplitude  to  the 
isotherm  perturbations,  and  produced  maximum  temperature  gradients  which 
are  approximately  twice  as  great  as  those  observed. 

As  was  the  case  with  the  surface  air  temperature,  the  largest  sim¬ 
ulation  errors  at  800  mb  (Fig.  17)  occur  over  the  continents,  where 


See  footnote,  p.  8  and  also  Fig.  18. 
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errors  of  5  to  10  deg  C  are  common  in  higher  latitudes.  Over  the  globe, 
the  simulated  average  800-mb  temperature  is  1  deg  C  below  the  observed 
January  value  of  4.1  deg  C.  The  large  simulation  errors  over  Antarctica 
(exceeding  26  deg  G)  are  similar  to  those  for  surface  air  temperature 
(Fig.  14),  inasmuch  as  the  surface  pressure  itself  is  below  800  mb. 

This  excessively  warm  air  is  due  to  the  models  use  of  an  erroneously 
low  albedo,  as  noted  above.  This  error  is  particularly  evident  in  the 
zonal  averages  of  Fig.  18,  where  the  anomalous  increase  of  temperature 
south  of  about  65°S  coincides  with  the  extent  of  the  ice  over  the  Ant¬ 
arctic  ocean  and  continent  (see  Figs.  2  and  3). 

Figure  20  records  the  departure  of  the  simulated  400-mb  temperature 
from  that  observed  and  reveals  (with  the  exception  of  Antarctica)  that 
the  largest  errors  occur  in  the  equatorial  and  subtropical  latitudes, 
in  contrast  to  the  error  patterns  at  both  the  surface  (Fig.  14)  and  at 
800  mb  (Fig.  17).  Over  the  Caribbean  and  over  the  Amazon,  for  example, 
the  simulated  40Q-mb  temperature  is  some  12  deg  C  above  that  observed 
for  January.  The  zonally  averaged  400-mb  temperature  in  Fig.  21  is 
about  6  deg  C  above  that  observed  between  40°N  and  20°S,  and  contri¬ 
butes  to  a  global  average  temperature  3.7  deg  above  the  observed  value 
of  24.2  deg  C.  Perhaps  of  more  significance  is  the  model’s  simulation 
of  a  midlatitude  meridional  temperature  gradient ,  which  is  about  60 
percent  too  great  in  both  hemispheres  (see  Fig.  21).  This  error  is 
reflected  in  the  mean  zonal  winds  discussed  below,  and  is  most  likely 
related  to  the  model’s  simulation  of  excessive  convective  precipitation 
in  the  lower  latitudes,  with  its  upper  level  receiving  too  much  (con¬ 
vective)  latent  heating. 

As  was  the  case  at  the  surface  and  at  800  mb,  the  air  temperature 
at  400  mb  over  Antarctica  is  systematically  higher  than  that  observed. 
South  of  about  70 °S  the  simulated  meridional  temperature  gradient  is 
of  the  wrong  sign,  with  a  10-deg  C  error  at  90°S  (Fig.  21).  This  dis¬ 
crepancy  is  due  to  the  albedo  error  noted  earlier. 

ZONAL  WIND  AT  800  MB  AND  400  MB 

The  simulated  January  zonal  (u)  wind  component  at  800  mb  and  400  mb 
is  shown. in  Figs.  22  and  23.  Because  of  the  lack  of  corresponding  ob¬ 
served  wind  data,  global  maps  of  the  observed  values  and  the  simulation 


Fig.  16b  —  The  observed  temperature  for  January  at  800  mb,  with  isolines  at  5-deg  C  intervals 
and  the  0-deg  C  isotherm  dashed.  Based  on  data  of  Crutcher  and  Meserve  (1970)  and 
Taljaard  et  al.  (1969)  as  summarized  by  Schutz  and  Gates  (1971), 
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Fig.  19a  —  The  simulated  temperature  for  January  at  400  mb,  with  isolines  every  5  deg  C  and 
the  -20-deg  C  isotherm  dashed. 
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The  observed  temperature  for  January  at  400  mb,  with  isolines  every  5  deg  C 
and  the  -20-deg  C  isotherm  dashed.  Based  on  data  of  Crutcher  and  Meserve  (1970) 
and  Taljaard  et  al.  (1969)  as  summarized  by  Schutz  and  Gates  (1971). 


Fig.  21  —  The  zonal  average  January  temperature  at  400  mb  as  simulated  (full  line)  and  as  observed 
(dashed  line),  based  on  the  data  of  Fig.  19. 
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errors  cannot  be  given.  In  general,  the  simulated  mean  zonal  winds 
clearly  show  the  midlatitude  jet  streams,  with  maximum  westerlies  off 
the  east  coasts  of  North  America  and  Asia,  and  with  the  westerly  winds 
generally  increasing  with  height.  These  simulations  correspond  reason¬ 
ably  well  with  the  observed  mean  geos trophic  zonal  wind,  which  may  be 
calculated  from  the  observed  geopotential  fields  (Figs.  7b  and  10b), 
and  with  the  available  observed  winds  in  the  northern  hemisphere  shown 
by  Schutz  and  Gates  (1971).  The  strongest  simulated  zonal  winds  near 
40°N  in  Figs.  22  and  23,  however,  are  approximately  twice  the  strength 
of  either  of  these  observed  measures .  In  the  lower  latitudes ,  the 
model  developed  weak  easterlies  between  about  25°N  and  20°S  at  800  mb, 
and  between  about  10°N  and  10°S  at  400  mb. 

In  spite  of  the  sparsity  of  the  observed  data,  a  useful  quantita¬ 
tive  comparison  can  be  made  with  zonal  averages  as  in  Figs.  24  and  25. 

At  800  mb  the  agreement  with  observation  may  be  considered  generally 
good,  although  the  southern  hemisphere  westerlies  and  the  northern 
hemisphere  easterlies  are  about  5  deg  too  far  north.  A  more  important 
simulation  error,  however,  is  the  systematic  overprediction  of  the  speed 
of  the  northern  hemisphere  westerlies;  the  simulated  800-mb  maximum 
zonal  wind  speed  of  14  m  sec  ^  at  40 °N  is  nearly  twice  the  8  m  sec 
speed  of  the  strongest  observed  mean  zonal  winds  (either  actual  or  geo- 
strophic) .  At  400  mb  the  speed  of  the  simulated  maximum  zonal  westerlies 
at  40°N  is  35  m  sec  \  a  67  percent  overprediction  of  Crutcher's  (1961) 
maximum  observed  zonal  wind  of  21  m  sec  ^  (at  35 °N)  and  a  25  percent 
overprediction  of  the  maximum  mean  zonal  geostrophic  wind  of  28  m  sec 
(Crutcher  and  Meserve,  1970).  The  maximum  westerlies  are  also  simulated 
between  5  and  10  deg  too  far  north.  These  errors  are  related  to  the 
excessively  strong  gradient  of  geopotential  noted  earlier,  which  in 
turn  is  related  to  the  excess  warmth  of  the  tropical  atmosphere.  The 
accompanying  increased  meridional  temperature  gradient  thereby  contri¬ 
butes  to  the  intensification  of  the  midlatitude  disturbances. 

It  is  also  of  interest  to  compare  the  simulated  and  observed  vert¬ 
ical  wind  shears.  From  Figs.  24  and  25,  the  maximum  difference  between 
the  simulated  zonal  winds  at  400  mb  and  800  mb  is  about  22  m  sec  at 
36°N,  and  about  19  m  sec"1  at  35°S.  The  maximum  shear  of  the  observed 
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zonal  winds  in  the  northern  hemisphere  between  400  mb  and  800  mb  occurs 
near  35°N  and  is  about  16  m  sec  ^  (Crutcher,  1961).  The  maximum  sim¬ 
ulated  shear  is  thus  about  37  percent  higher  than  that  observed  (in 
the  northern  hemisphere) .  In  terms  of  the  vertical  shear  of  the  ob¬ 
served  geostvophia  zonal  winds,  the  maximum  differences  between  400  mb 
and  800  mb  are  about  20  m  sec  ^  (at  36°N)  and  about  13  m  sec  ^  (at  48°S), 
according  to  the  data  of  Crutcher  and  Meserve  (1970)  as  summarized  by 
Schutz  and  Gates  (1971).  These  agree  quite  well  with  the  model’s  cal¬ 
culations  which  may  therefore  be  said  to  correctly  simulate  the  ob¬ 
served  thermal  wind  balance  in  the  troposphere. 

The  most  direct  way  to  reduce  the  excessive  strength  of  the  simu¬ 
lated  zonal  winds  (while  possibly  preserving  the  model’s  nearly  correct 
vertical  wind  shear)  would  seem  to  be  either  to  increase  the  surface 
friction  (see  Gates,  1972)  or  to  introduce  lateral  eddy  viscosity  that 
is  not  now  used  in  the  model  (see  Gates  et  al.,  1971). 

MERIDIONAL  WIND  AT  800  MB  AND  400  MB 

The  simulated  January  meridional  (v)  wind  component  at  800  mb  and 
at  400  mb  is  shown  in  Figs.  26  and  27.  As  with  the  zonal  wind,  avail¬ 
able  data  are  not  sufficient  to  permit  a  detailed  global  comparison  with 
observation.  In  general,  the  mean  meridional  winds  have  maximum 
strength  in  connection  with  the  standing  long  waves  in  the  middle  lati¬ 
tudes  of  the  northern  hemisphere.  Here  the  simulated  winds  agree 
reasonably  well,  in  both  position  and  magnitude,  with  the  available 
observations  in  the  northern  hemisphere  and  with  the  global  geos trophic 
winds  as  given  by  Schutz  and  Gates  (1971).  Generally,  the  pattern  of 
the  mean  meridional  wind  shows  coherence  between  the  levels  400  mb  and 
800  mb,  with  approximately  twice  as  great  speed  at  the  upper  level. 

The  prominent  feature  is  the  one  associated  with  the  Icelandic  low, 
where  the  mean  January  meridional  wind  exceeds  20  m  sec  at  400  mb. 

The  mean  meridional  winds  in  the  tropics  and  in  most  of  the  southern 
hemisphere  are  relatively  weak,  although  there  is  still  a  noticeable 
resemblance  between  the  winds  at  both  levels. 

The  zonal  averages  of  the  meridional  winds  present  a  somewhat  dif¬ 
ferent  picture  (Figs.  28  and  29).  Here  the  simulated  mean  meridional 


The  zonal  wind  at  800  mb  simulated  for  January.  The  isolines  are  at  intervals  of 
5  m  sec“l,  with  the  zero  line  dashed.  The  symbols  W  and  E  denote  local  maxima  of 
winds  from  the  west  and  east. 


Latitude  (deg) 


•  24  The  zonal  average  of  the  mean  zonal  winds  at  800  mb  simulated  for  January  (full  line) , 
together  with  the  distribution  of  the  corresponding  observed  geos trophic  winds  (dashed 
line)  as  given  by  Schutz  and  Gates  (1971)  based  on  data  of  Crutcher  and  Meserve  (1970) 
and  Taljaard  et  al»  (1969).  Also  shown  by  the  open  circles  in  the  northern  hemisphere 
are  the  zonal  averages  of  the  observed  zonal  wind  (for  December,  January,  February) 
from  Schutz  and  Gates  (1971)  based  on  data  of  Crutcher  (1961). 
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Fig.  25  —  The  zonal  average  of  the  mean  zonal  winds  at  400  mb  simulated  for  January  (full  line), 
together  with  the  distribution  of  the  corresponding  observed  geostrophic  winds. (dashed 
line)  as  given  by  Schutz  and  Gates  (1971)  based  on  data  of  Crutcher  and  Meserve  (1970) 
and  Taljaard  et  al.  (1969).  Also  shown  by  the  open  circles  in  the  northern  hemisphere 
are  the  zonal  averages  of  the  observed  zonal  wind  (for  December,  January,  February) 
from  Schutz  and  Gates  (1971)  based  on  data  of  Crutcher  (1961). 


Meridional  wind  at  800  mb  (m  sec”1) 


Latitude  (deg) 

Fig.  28  —  The  zonal  average  of  the  mean  meridional  winds  at  800  mb  simulated  for  January 

(full  line) ,  with  the  corresponding  observed  winds  (dashed  line)  in  the  northern 
hemisphere.  From  Schutz  and  Gates  (1971),  based  on  data  of  Crutcher  (1961). 


Meridional  wind  at  400  mb  (m  sec”1) 


1.5 


Fig.  29  —  The  zonal  average  of  the  mean  meridional  winds  at  400  mb  simulated  for  January 

(full  line),  with  the  corresponding  observed  winds  (dashed  line)  in  the  northern 
hemisphere.  From  Schutz  and  Gates  (1971),  based  on  data  of  Crutcher  (1961). 
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flow  is  dominated  by  the  components  of  the  Hadley  circulation  in  the 
tropics,  where  the  zonally  averaged  meridional  wind  exceeds  1  m  sec 
A  maximum  southward  motion  is  simulated  between  20 °N  and  30°N  at  800  mb 
(Fig.  28),  and  a  nearly  equal  maximum  northward  flow  at  400  mb  (Fig.  29) 
at  the  same  latitudes.  A  reversed  and  slightly  weaker  circulation  is 
simulated  in  the  southern  hemisphere  between  20°S  and  30°S  at  both 
levels.  These  motions  are  essentially  nongeos trophic ,  since  the  zonally 
averaged  meridional  geostrophic  wind  vanishes  identically.  In  comparison 
with  the  zonal  mean  of  the  available  observed  meridional  winds  for  the 
northern  hemisphere  [as  given  by  Schutz  and  Gates  (1971),  based  on  data 
of  Crutcher  (1961)],  the  simulated  circulation  is  systematically  dis¬ 
placed  toward  the  north.  At  800  mb  (Fig.  28)  the  maximum  equatorward 
flow  is  simulated  some  20  deg  north  of  the  observed  maximum  at  5°N,  while 
at  400  mb  the  maximum  poleward  flow  is  simulated  about  10  deg  too  far 
north.  A  northward  displacement  is  also  evident  in  other  features  of 
the  mean  meridional  winds,  although  the  simulated  speeds  are  in  reason¬ 
ably  good  agreement  with  those  observed  (in  the  northern  hemisphere) . 

VERTICAL  M)TI0N  AND  MASS  FLUX 

The  mean  meridional  velocities  at  400  mb  and  800  mb,  however,  do 
not  portray  the  mass  flux  associated  with  the  mean  meridional  circula¬ 
tion  in  the  model,  because  of  the  varying  pressure  at  the  earth’s  sur¬ 
face.  To  examine  this  flux,  we  may  write  the  continuity  equation  in 
the  form 


9t r 
3t 


£  V 


7TV  +  (to)  *  0 


(7) 


where  the  symbols  are  as  previously  defined.  Denoting  a  time  average 
by  (  )  and  a  zonal  average  by  [  ] ,  this  becomes 
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assuming  (9[7r]/9t)  =  0  over  the  (30-day)  time  period.  A  streamfunction 
for  the  total  mass  flux  may  then  be  defined  for  which 
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2Tr'a^cos  cp  [ tto 3 / g  =  -3i|>/3cp 

(9a) 

27r'a  cos  cp  [ttv] /g  =  dty/do 

(9b) 

where  a  is  the  earth1  s  radius,  cp  is  latitude,  g  is  gravity,  and  tt*  =  3.14... 
(to  avoid  confusion  with  tt  =  pg  -  p^  previously  defined)  . 

We  may  use  either  of  these  expressions  to  find  ip,  depending  upon 
the  data  available.  As  suggested  by  Mints,  I  have  used  (9a)  and  found 
Jp  by  an  integration  with  respect  to  latitude.  Between  two  latitudes 
and  Cp^,  we  thus  have 


cos  cp  [tt a2]  dcp 


(10) 


With  the  condition  ip  *  0  at  the  poles  and  at  both  O  *  0  and  C*l,  this 
expression  was  used  to  calculate  the  streamf unction  distribution  shown 
in  Fig.  30b.  The  values  of  I used  for  this  purpose  are  shown  in 
Fig.  30a,  as  given  by  the  model’s  mean  January  simulation.  Since  ip 
is  found  only  at  the  model’s  midlevel,  the  simulated  streamf unction  is 
symmetric  about  0=1/2.  A  comparison  with  the  observed  mass  flux 
streamf unction  given  by  Newell  et  al.  (1970),  also  shown  in  Fig.  30b, 
reveals  that  the  simulated  mean  meridional  circulation  has  a  number  of 
realistic  features. 

We  note  in  particular  that  the  strength  of  the  simulated  northern- 

12  —1 

hemisphere  Hadley  circulation  (152  x  10  g  sec  )  is  close  to  that  given 
by  Newell  et  al.  (1970)  for  December,  January,  and  February,  although 
it  is  somewhat  below  the  values  180  x  10^  g  sec  ^  and  230  x  10  g  sec 
given  by  Oort  and  Rasmusson  (1970)  for  January  and  by  Palmen  and  Newton 
(1969)  for  the  northern  winter.  The  position  of  the  simulated  northern 
Hadley  cell,  however,  is  approximately  15  deg  too  far  north.  The  model’s 
northern  Ferrel  cell  is  too  weak  and  is  much  narrower  than  that  observed, 
although  it  is  simulated  at  approximately  the  correct  mean  latitude. 

This  may  be  due  to  the  model’s  tendency  to  underestimate  the  strength 


Vertical  motion  Tr&2  ,  (mb  hr-1) 
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Fig.  30a  —  The  zonal  average  of  the  mean  vertical  motion  simulated  for  January,  as  given 
by  the  measure  ttc 2 .  Neglecting  surface  pressure  variations,  1  mb  hr~l  here 
corresponds  to  a  vertical  velocity  of  about  0.4  cm  sec~^  at  600  mb.  The 
shading  denotes  regions  of  upward  vertical  motion  (see  also  Fig.  30b). 


Pressure  (mb)  Sigma  (a) 


Fig.  30b  The  streamf unction  for  the  time-mean,  zonally  averaged  mass  flux  in  10^  g  sec-l# 
That  simulated  for  January  is  shown  above  (symmetric  about  a  =  1/2),  and  that 
observed  (for  December,  January,  February)  is  shown  below  as  given  by  Newell  et  al. 
(1970).  The  direction  of  flow  is  given  by  the  arrows. 
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of  the  transient  midlatitude  cyclones  noted  earlier,  with  which  the 
bulk  of  the  nongeos trophic  transport  in  the  Ferrel  cell  may  be 
associated. 

In  the  southern  hemisphere,  the  simulated  strength  (and  width)  of 
the  Hadley  circulation  is  much  greater  than  that  observed,  at  least 
according  to  Newell  et  al.  (1970).  The  cause  of  this  discrepancy  is 
not  clear,  although  it  is  interesting  to  note  that  the  southern  Ferrel 
cell  (as  well  as  the  weaker  cells  south  of  it)  is  also  simulated  with 
too  great  a  strength.  It  may  be  that  the  model  systematically  over¬ 
estimates  the  strength  of  the  southern  hemisphere  circulation  for  some 
reason  (such  as  the  use  of  too  low  a  sea-surface  temperature) ,  or  that 
the  available  observations  in  the  southern  hemisphere  are  insufficient 
to  determine  the  mean  meridional  mass  flux  adequately* 

After  the  mass  flux  has  been  determined,  the  mean  mass-weighted 
meridional  velocity  between  0=1  and  0=1/2  may  be  found  from  (9b)  as 

_  I  ™ —  i 

[v  ]  =* — — —  =  -2iJ;g(2'fr,a  cos  cp  [tt] )  (11) 

J  [tt] 

with  lv±]  between  the  levels  0=1/2  and  o  =  0  (or  200  mb)  given  by  a 
change  of  sign.  In  the  northern  Hadley  cell,  for  example,  we  thus  find 
a  maximum  [v^]  =*  — 1.0  m  sec  ^  at  about  20°N  in  contrast  to  the  maximum 
value  of  about  -1.3  m  sec  ^  found  at  30°N  at  800  mb  (Fig.  28).  The 
correlation  between  IT  and  v  thus  not  only  insures  [ttv^]  +  [ttv^]  = 
but  also  permits  an  effective  mass-flux  velocity  different  from  the 


mean  meridional  velocities  themselves . 
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V.  THE  SIMULATED  AVERAGE  MOISTURE  AND  PRECIPITATION 


In  addition  to  the  fields  of  the  primary  variables  of  pressure, 
temperature,  and  wind,  the  distribution  of  the  atmospheric  water  vapor 
and  of  its  sources  and  sinks  (evaporation  and  precipitation)  is  an 
important  measure  of  the  models  ability  to  simulate  climate.  In  this 
section  the  global  distributions  of  these  variables  as  simulated  for 
January  in  the  control  integration  will  be  considered  and  compared  with 
the  available  observations. 

RELATIVE  HUMIDITY  AT  800  MB 

The  simulated  January  relative  humidity  at  800  mb  is  given  in 
Fig.  31,  along  with  the  corresponding  observed  pattern.  Here  the  sim¬ 
ulated  relative  humidity  was  determined  from  the  computed  mixing  ratio  4 
and  the  temperature  T^  at  the  level  3  as  described  in  Gates  et  al . 
(1971) ,  whereas  the  observed  data  are  those  given  by.  Schutz  and  Gates 
(1971)  based  on  the  dew-point  data  of  Crutcher  and  Meserve  (1970)  and 
Taljaard  et  al.  (1969). 

Only  in  a  gross  sense  can  the  simulated  field  (Fig.  31a)  be  con¬ 
sidered  a  satisfactory  approximation  to  the  observed  humidity  field 
(Fig.  31b).  The  prominent  observed  features  are  the  dry  regions  over 
Africa,  India,  and  Australia,  and  in  the  southwestern  portions  of  North 
America,  South  America,  and  South  Africa.  The  model  has  reproduced  the 
low  relative  humidity  over  these  areas  with  reasonable  accuracy.  The 
regions  of  higher  humidity,  however,  are  not  well  depicted  in  general, 
and  there  are  many  dissimilar  features  in  the  simulated  and  observed 
patterns.  Particularly  noticeable  is  the  model’s  failure  to  simulate 
the  regions  of  high  humidity  over  South  America,  Africa,  and  Indonesia. 

The  errors  of  the  simulated  800-mb  relative  humidity  shown  in 
Fig.  32  are  generally  largest  in  the  tropics,  where  local  errors  of  40 
percent  are  not  uncommon.  We  note,  moreover,  that  here  the  predominant 
error  is  negative;  i.e.,  the  simulation  is  drier  than  observation.  Part 
of  this  error  may  be  due  to  the  model’s  use  of  a  fictitious  wind  at 
a  *  7/8  to  accomplish  the  horizontal  moisture  advection  in  (4) ,  and  part 
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Fig  31b  —  The  observed  relative  humidity  for  January  at  800  mb,  with  isolines  at  10-percent 
intervals  and  the  60-percent  isoline  dashed.  From  Schutz  and  Gates  (1971) ,  based 
on  data  of  Crutcher  and  Meserve  (1970)  and  Taljaard  et  al.  (1969). 
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Fig.  32  —  The  difference  between  the  simulated  and  observed  800— mb  relative  humidity  (Fig.  31a 
minus  Fig.  31b).  The  isoline  interval  is  10  percent  with  the  zero  line  dashed. 


Relative  humidity  at  800  mb  (%) 
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Fig.  33  —  The  zonal  average  of  the  mean  800  mb  relative  humidity  simulated  for  January 
(full  line) ,  with  the  corresponding  observed  distribution  (dashed  line) , 
based  on  the  data  of  Fig.  31. 


-62- 


may  be  due  to  the  assumption  of  a  completely  dry  atmosphere  above  the 
midlevel  (approximately  600  mb) ,  and  the  consequent  absence  of  vertical 
moisture  advection.  Systematic  errors  in  the  simulated  moisture- 
addition  rate  (E  -  C)  may  also  contribute,  as  will  be  discussed  below. 

This  low-latitude  dryness  is  clearly  seen  in  terms  of  the  zonal 
averages  given  in  Fig.  33.  Poleward  of  about  30°N  and  30°S,  however, 
the  simulation  has  substantially  overestimated  the  relative  humidity, 
resulting  in  the  model's  global  average  relative  humidity  of  59  per¬ 
cent  approaching  the  observed  value  of  58  percent.  This  suggests  that 
the  model  may  be  so  efficient  in  the  poleward  advection  of  moisture 
that  the  higher  latitudes  are  oversupplied  while  the  tropics  are  de¬ 
pleted.  As  is  discussed  below,  it  is  also  possible  that  the  model 
gives  too  little  large-scale  precipitation  in  the  middle  latitudes 
(and  thereby  does  not  deplete  the  midlatitude  atmospheric  moisture 
supply  sufficiently) ,  whereas  it  gives  too  much  convective  precipita¬ 
tion  in  the  lower  latitudes  (and  thereby  overdepletes  the  tropical 
atmospheric  moisture) . 

CLOUDINESS 

Clouds  in  the  Mintz-Arakawa  two-level  model  are  classified  into 
three  types,  depending  upon  the  particular  condensation  processes  in¬ 
volved.  Cloud  type  1  (CL1  in  the  model)  represents  deep  convective 
cloudiness  at  level  1  and  is  assumed  to  accompany  either  midlevel  or 
penetrating  convection  (which,  in  turn,  depends  upon  the  vertical 
stability) .  A  second  form  of  convective  cloud  simulated  is  type  3 
(CL3  in  the  model) ,  which  represents  shallow  convective  cloudiness  or 
cumulus  at  level  3.  Both  types  of  cloudiness  are  determined  from  the 
empirical  formula 


CL1 ,  CL3  =  -1.3  +  2.6  RH3  (12) 

where  RH^  is  the  relative  humidity  at  level  3  (as  in  Fig.  31a),  with 
the  additional  stipulation  that  0  £  (CL1,  CL3)  ^  1.  These  convective 
cloud  forms  occur  only  when  the  model's  criteria  for  convective  insta¬ 
bility  are  met,  after  which  there  is  an  adjustment  of  both  temperature 
and  humidity  to  restore  stability. 
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The  third  cloud  form  simulated  is  the  large-scale  cloudiness  asso¬ 
ciated  with  nonconvective  large-scale  condensation.  In  the  model  this 
cloudiness  (CL2)  is  set  equal  to  1  if  large-scale  precipitation  is 
occurring,  which  itself  depends  on  the  apparent  supersaturation  at  level 
3  as  a  result  of  vertical  motion.  If  large-scale  precipitation  (or  con¬ 
densation)  is  not  occurring,  then  CL2  =  0.  Thus,  middle  cloudiness  is 
here  simulated  as  either  overcast  or  clear  sky. 

Each  of  the  above  cloud  types  affects  the  simulated  absorption  of 
both  short-  and  long-wave  radiation  in  the  model  atmosphere  as  described 
elsewhere  (Gates  et  al.,  1971).  The  fraction  of  the  sky  assumed  to  be 
cloud-covered  is  given  by  the  expression  CL  =  max  (CL1,  CL2,  CL3) ,  and 
it  is  this  total  cloudiness  measure  whose  simulated  global  distribution 
is  shown  in  Fig.  34.  The  observed  January  cloudiness  shown  in  Fig.  35a 
is  that  given  by  Schutz  and  Gates  (1971)  for  the  northern  hemisphere, 
as  prepared  from  data  furnished  by  the  Environmental  Technical  Appli¬ 
cations  Center  of  the  U.S.  Air  Force  (1971). 

In  general,  the  January  simulation  has  specified  the  maximum  cloud- 
iess  over  the  North  Pacific  and  North  Atlantic  oceans  with  reasonable 
accuracy.  The  model’s  possible  ability  to  depict  the  relatively  clear 
skies  over  the  southwestern  United  States,  over  North  Africa,  and  over 
eastern  Asia  is  obscured  by  its  apparent  tendency  to  specify  a  low 
amount  of  cloudiness  over  nearly  all  continental  areas.  This  same 
tendency  is  responsible  for  the  simulation  of  relatively  low  cloudiness 
amounts  over  Europe,  where  there  is  substantial  cloudiness  observed 
(Fig.  35a). 

These  errors  have  resulted  in  differences  between  the  simulated 
and  observed  cloudiness  of  as  much  as  4/10  over  widespread  areas,  as 
shown  in  Fig.  35b.  In  terms  of  zonal  averages,  shown  in  Fig.  36,  the 
model’s  mean  cloudiness  is  systematically  lower  than  the  northern- 
hemisphere  observations  by  approximately  2/10,  although  the  cloudiness 
maximum  is  reproduced  in  about  the  correct  midlatitude  position.  Also 
shown  in  Fig.  36  are  the  zonal  averages  of  the  TIROS  satellite  cloudiness 
observations  of  Miller  (1970)  ,  here  taken  from  the  compilation  of  Schutz 
and  Gates  (1971).  Even  though  the  satellite  data  contain  systematic 
errors  due  to  the  lack  of  discrimination  between  clouds  and  highly- 
reflective  ground  surfaces  such  as  snow,  ice,  and  desert,  the  model’s 
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Fig.  34  —  The  simulated  total  cloudiness  for  January,  with  isolines  at  0.1-cloudiness  intervals 
and  the  0.5  isoline  dashed. 
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35a  —  The  observed  January  cloudiness  in  the  northern  hemisphere,  with  isolines  at  0.1- 
cloudiness  intervals  and  the  0.5  isoline  dashed.  From  Schutz  and  Gates  (1971), 
as  prepared  from  data  supplied  by  the  Environmental  Technical  Applications  Center 
of  the  U.S.  Air  Force  (1971). 


35b  —  The  difference  between  the  simulated  and  observed  cloudiness  (Fig.  34  minus  Fig.  35a) 
in  the  northern  hemisphere,  with  the  isolines  at  0.2-cloudiness  intervals  and  the 
zero  isoline  dashed. 


Total  cloudiness  (fractions) 


Latitude  (deg) 

Fig.  36  —  The  zonal  average  of  the  mean  cloudiness  as  simulated  for  January  (full  line), 
together  with  the  observed  distribution  (dashed  line).  From  Schutz  and  Gates 
(1971),  based  on  satellite  data  of  Miller  (1970).  Also  shown  in  the  northern 
hemisphere  (open  circles)  are  the  observed  values  (from  both  satellite  and 
conventional  observations)  based  on  the  data  of  Fig.  35.  The  symbol  L  de¬ 
notes  the  simulated  low  cloudiness  only. 
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January  cloudiness  is  clearly  too  low,  with  a  simulated  global  average 
of  0,26  compared  to  an  observed  value  of  0.41  from  Miller’s  data. 

Most  of  the  model’s  cloudiness  is  here  due  to  the  simulation  of 
low  clouds  (type  3  or  CL3  described  above) ,  which  are  associated  pri¬ 
marily  with  convective  heating  over  the  oceans;  the  zonal  average  cloud¬ 
iness  due  to  this  cloud  type  alone  is  also  shown  in  Fig.  36.  Only  north 
of  about  50 °N  does  the  cloudiness  associated  with  large-scale  precipita¬ 
tion  contribute  as  much  as  2/10  cloudiness  to  the  zonal  means.  At  no 
latitude  does  the  simulated  high  cloudiness  (CL1)  give  a  zonal  average 
greater  than  5/10.  This  suggests  that  it  may  be  desirable  to  increase 
systematically  the  total  cloudiness  simulated  by  the  model  at  all  lati¬ 
tudes  ,  and  in  particular  to  increase  the  middle-level  cloudiness  whether 
or  not  it  is  associated  with  precipitation. 

PRECIPITATION 

The  precipitation  rate  in  the  moisture-balance  equation  (4)  is 
taken  equal  to  the  total  rate  of  condensation  due  to  both  large-scale 
and  convective  processes.  The  precipitation  associated  with  large- 
scale  motions  is  given  by  the  apparent  supersaturation  (due  primarily 
to  large-scale  vertical  motion) ,  while  the  convective  precipitation  is 
equal  to  the  condensation  resulting  from  the  parameterized  midlevel 
and  penetrating  convection  (due  to  vertical  instability).  The  total 
precipitation  rate  simulated  for  January  is  shown  in  Fig.  37,  together 
with  the  observed  distribution  for  December,  January,  and  February  as 
given  by  Schutz  and  Gates  (1972)  from  data  of  Moller  (1951) .  Here  the 
average  simulated  precipitation  rate  (Fig.  37a)  was  determined  from 

* 

the  solutions  every  1/2  hour  (the  minimum  sampling  interval  possible) , 
rather  than  from  the  solutions  every  6  hours.  Sampling  less  frequently 
produces  a  distorted  rainfall  distribution  in  the  tropics  where  con¬ 
vective  precipitation  is  predominant  (Gates,  1972). 

A  comparison  of  Figs.  37a  and  37b  shows  that  the  model  has  simu¬ 
lated  the  arid  regions  southwest  of  the  North  American,  South  American, 

The  simulated  precipitation  distribution  (Fig.  37a)  is  from  a 
separate  integration  (see  Gates,  1972),  rather  than  from  the  control 
integration  used  elsewhere  in  this  report. 


Fig.  37a  —  The  simulated  precipitation  rate  for  January  (in  a  modified  version  of  the  control 
integration)  ,  with  isolines  at  intervals  of  1  mm  day-1  and  the  2  mm  day-1  isoline 
dashed . 
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Fig.  37b  —  The  observed  precipitation  for  December,  January,  and  February,  with  isolines  at 
intervals  of  1  mm  day“^  and  the  2  mm  day“l  isoline  dashed.  Based  on  data  of 
Moller  (1951)  as  summarized  by  Schutz  and  Gates  (1972). 


38  —  The  difference  between  the  simulated  and  observed  precipitation  rate  (Fig.  37a 

minus  Fig.  37b).  The  isoline  interval  is  1  mm  day-1  with  the  zero  isoline  dashed 


Fig.  39  —  The  zonal  average  of  the  precipitation  rate  as  simulated  for  January  (full  line) 
and  as  observed  (dashed  line),  based  on  the  data  of  Fig.  37.  The  dash-dotted 
line  gives  the  observed  precipitation  from  the  data  of  Lvovitch  and  Ovtchinnikov 
(1964)  as  modified  by  Schutz  and  Gates  (1971)  for  December,  January,  and  February. 
The  symbol  L  denotes  the  zonal  averages  of  the  simulated  large-scale  precipitation 
rate  only,  and  the  open  circles  denote  that  portion  of  the  total  simulated  pre¬ 
cipitation  that  falls  over  the  oceans. 
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African,  and  Australian  continents  with  reasonable  fidelity.  The  most 
serious  errors  occur  in  the  specification  of  the  areas  of  maximum  pre¬ 
cipitation  over  the  North  Atlantic  (where  the  simulated  average  pre¬ 
cipitation  rate  of  15  mm  day  ^  is  several  times  that  observed) ,  and  in 
the  tropics  (where  the  simulated  rainfall  rate  is  generally  about  twice 
that  observed  and  displaced  poleward) .  We  may  also  note  the  near 
absence  of  simulated  precipitation  over  the  bare  land  areas  (with  the 
notable  exceptions  of  South  America  and  Africa) .  This  feature  results 
in  a  poor  representation  of  the  observed  precipitation  over  eastern 
North  America,  Europe,  and  (especially)  over  Southeast  Asia.  Although 
there  is  some  uncertainty  in  Moller’s  (1951)  observed  precipitation 
distribution,  the  global  rainfall  given  by  Lvovitch  and  Ovtchinnikov 
(1964)  (see  Schutz  and  Gates,  1971)  shows  a  similar  general  distribution. 
The  precipitation  simulation  errors  shown  in  Fig.  38  may,  therefore,  be 
considered  representative  of  the  model. 

In  the  zonal  averages  shown  in  Fig.  39,  the  prominent  error  is 
again  the  model’s  excessive  precipitation  between  approximately  25 °N 
and  25°S.  The  bulk  of  this  precipitation  falls  on  the  oceans,  with 
a  maximum  at  18°N  and  the  subtropical  dry  zone  in  the  northern  (winter) 
hemisphere  displaced  some  15  deg  northward  of  its  observed  position. 

The  simulated  higher-latitude  precipitation  is  also  too  great  (in  the 
northern  hemisphere) ,  although  the  maximum  is  correctly  positioned  at 
about  40°N.  This  error  is  also  partly  due  to  the  simulation  of  exces¬ 
sive  convection,  as  large-scale  and  convective  condensation  contribute 
approximately  equally  to  the  total  precipitation  poleward  of  30°N  (see 
Fig.  39).  [The  separate  global  patterns  of  the  convective  and  large- 
scale  precipitation  simulated  in  this  experiment  have  been  given  by 
Gates  (1972).]  Only  in  the  southern  hemisphere  between  50°S  and  65°S 
is  the  precipitation  significantly  underestimated  by  the  model,  although 
here  the  observations  themselves  may  be  less  reliable. 

Another  interesting  feature  of  Figs.  37  and  39  is  the  relatively 
heavy  precipitation  simulated  over  Antarctica.  In  contrast  to  the 
generally  stable  and  dry  atmosphere  observed  between  70°S  and  90°S, 
the  model  produced  appreciable  deep  convection  in  this  region.  This 
is  another  result  of  the  inadvertent  surface-albedo  error  previously 
noted. 
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On  a  global  average,  the  simulated  total  precipitation  rate  of 
4,1  mm  day  ^  (of  which  3.5  mm  day  ^  is  convective)  is  some  86  percent 
higher  than  the  observed  rate  of  2.2  mm  day  from  Moller's  data  in 
Fig.  39.  As  is  noted  earlier,  nearly  all  of  this  error  is  due  to  the 
excessive  convective  precipitation  simulated  in  the  tropics.  The  ex¬ 
cessive  latent  heat  released  by  this  parameterized  convection  may  also 
contribute  to  the  model's  specification  of  too  high  a  temperature  at 
400  mb  previously  cited  (see  Fig.  21),  as  well  as  to  the  low  tropical 
relative  humidity  and  cloudiness  (see  Figs.  33  and  36). 


EVAPORATION 

The  evaporation  rate,  E,  in  the  moisture-continuity  equation  (4) 
is  determined  from  the  expression 


E  =  P^WCqg  -  q4)pw 


(13) 


where  is  the  surface  (level  4)  air  density,  W  is  a  measure  of  the 
surface  wind  speed,  is  the  density  of  water,  is  the  surface  drag 
coefficient  (itself  dependent  upon  the  surface  wind  speed) ,  and  and 
q^  are  the  mixing  ratios  at  the  ground  and  in  the  surface  air  (just 


above  the  ground) . 


Over  the  oceans  q  is  taken  as 
g 


appropriate  to  the  assigned  sea-surface  temperature 


the  saturation  value 
(see  Fig.  3),  while 


over  land  (and  ice)  surfaces  q  is  determined  from  the  local  moisture 


and  heat  balance.  Over  all  surfaces,  q^  is  principally  determined  by 
the  moisture  at  level  3  (q^) ,  although  it  too  is  affected  by  the  surface 
heat  and  moisture  balances.  For  details  of  these  processes,  the  reader 


is  referred  to  Gates  et  al.  (1971). 

The  average  evaporation  rate  simulated  for  January  is  shown  in 
Fig.  40,  together  with  the  observed  distribution  from  Schutz  and  Gates 
(1971)  based  on  data  of  Budyko  (1963).  We  see  at  once  that  the  overall 
pattern  of  the  observed  evaporation  rate  is  simulated  reasonably  well 
by  the  model,  although  the  rate  is  overestimated  by  as  much  as  100  per¬ 
cent  in  the  tropical  latitudes.  Largely  through  the  controlling  in¬ 
fluence  of  the  sea-surface  temperature,  the  model  has  simulated  the 
zone  of  lower  evaporation  near  the  equator  in  the  Indian  and  Pacific 
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oceans,  but  has  even  here  overestimated  the  rate  in  absolute  terms. 

Over  land,  relatively  little  evaporation  is  simulated  except  where  the 
groundwater  level  is  kept  high  by  heavy  rainfall.  In  the  model  this 
occurs  in  January  only  over  the  southern  Amazon  Basin,  over  south- 
central  Africa,  and  in  some  coastal  regions.  In  the  difference  map 
shown  as  Fig.  41,  the  largest  simulation  errors  are  found  off  the  east 
coast  of  the  United  States  and  off  the  west  coast  of  Central  America. 

Here  the  computed  average  evaporation  rates  are  several  times  those 
observed,  with  the  error  occasionally  exceeding  10  mm  day 

This  apparent  tendency  of  the  model  to  overpredict  the  evaporation 
is  perhaps  more  clearly  seen  in  the  zonal  averages  shown  in  Fig.  42. 
Between  about  30°N  and  30 °S  the  general  pattern  of  the  simulated  evapor¬ 
ation  is  similar  to  that  observed,  but  the  computed  rates  are  approxi¬ 
mately  50  percent  greater  than  the  observed  ones.  Although  no  observa¬ 
tions  are  available  south  of  60°S,  the  high  rate  of  evaporation  simulated 
over  Antarctica  may  be  attributed  to  the  surface  albedo  error  discussed 
earlier. 

Averaged  globally,  the  simulated  January  evaporation  rate  of 
4.0  mm  day  ^  is  consistent  with  the  simulated  global  average  precipita¬ 
tion  rate  of  4.1  mm  day  ^  (Figs.  37  and  39).  Both  rates,  however, 
substantially  exceed  the  observed  global  average  rates  of  evaporation 
and  precipitation  (3.1  mm  day  ^  and  2.2  mm  day  .  The  model’s  over¬ 
prediction  of  the  precipitation  is  therefore  passed  on  in  the  form  of 
excessive  (and  approximately  equal)  evaporation.  It  is  apparent, 
however,  that  the  ’’observed”  long-term  average  global  rates  of  evapor¬ 
ation  and  precipitation  do  not  balance  as  well  as  do  the  simulated  ones. 
It  would  appear  that  the  mean  December,  January,  February  global  pre¬ 
cipitation  is  seriously  underestimated  by  Moller  (1951) ,  or  that  the 
evaporation  is  systematically  overestimated  in  the  data  of  Budyko 
(1963),  or  both.  Since  the  data  of  Lvovitch  and  Ovtchinnikov  (1964), 
as  summarized  by  Schutz  and  Gates  (1971) ,  give  a  global  mean  precipi¬ 
tation  rate  of  2.9  mm  day  \  the  greater  observational  error  would 
seem  to  lie  in  Holler’s  data. 

The  details  of  the  model’s  simulation  of  the  net  rate  of  moisture 
addition,  E  -  C,  are  given  elsewhere  (Gates,  1972),  Here  we  note. 
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however,  that  this  simulated  difference  is  positive  only  between  20° 
and  40°  latitude  in  each  hemisphere,  which  agrees  reasonably  well  with 
the  observed  average  moisture  source  for  January  from  Figs.  39  and  42. 

GROUND  WETNESS 

To  complete  this  presentation  of  the  model’s  hydrologic  cycle,  the 
mean  simulated  January  ground  wetness  is  shown  in  Fig.  43.  As  described 
in  detail  elsewhere  (Gates  et  al.,  1971),  the  ground  wetness  is  a 
dimensionless  measure  of  the  surface  wetness,  scaled  from  0  to  1.  Over 
water,  snow,  and  ice,  it  is  taken  as  unity,  whereas  over  land  it  is 
calculated  from  the  surface  water  balance.  When  the  (local)  precipita¬ 
tion  exceeds  evaporation,  the  ground  wetness  is  increased  toward  unity, 
with  a  portion  of  the  excess  water  assumed  lost  in  surface  runoff. 

When  evaporation  exceeds  precipitation,  the  local  ground  wetness  is 
decreased  (without  runoff). 

Although  there  are  no  observations  of  average  ground  wetness 
readily  available,  the  areas  of  dry  land  surface  in  the  simulated 
January  ground  wetness  shown  in  Fig.  43  correspond  roughly  with  the 
locations  of  the  world's  deserts.  The  most  notable  exceptions  are 
central  Africa,  India,  and  Indochina,  over  which  there  is  little  sim¬ 
ulated  rainfall  (Fig.  37a).  The  ground  wetness  itself  is  used  in  the 
calculation  of  the  ground  temperature  and  the  evaporation,  and  the 
dry  areas  in  Fig.  43  are  all  regions  of  low  simulated  evaporation 
(Fig.  40a).  The  only  areas  of  extensive  wet  ground  and  appreciable 
runoff  are  in  the  Amazon  Basin,  in  southeast  Africa,  and  in  Indonesia, 
which  correspond  well  with  the  regions  of  high  precipitation  observed 
in  January  as  shown  in  Fig.  37b. 


Fig.  42  —  The  zonal  average  of  the  evaporation  rate  as  simulated  for  January  (full  line) 
and  as  observed  (dashed  line),  based  on  the  data  of  Fig.  40. 


The  simulated  ground  wetness  in  January,  with  isolines  every  0.2  and  the  0.6 
isoline  dashed.  The  oceans  and  all  ice-covered  surfaces  are  considered  completely 
wet  (ground  wetness  El). 


-81- 


VI.  THE  SIMULATED  AVERAGE  SURFACE  HEATING 


Of  the  many  processes  which  contribute  to  the  net  heating  rate,  H, 
in  the  model,  only  a  few  can  be  readily  compared  at  present  with  corre¬ 
sponding  observed  global  distributions.  Here  we  shall  consider  only 
the  flux  of  sensible  heat  and  the  net  surface  absorption  of  shortwave 
radiation,  both  of  which  are  components  of  the  surface  heat  balance. 

The  heating  at  tropospheric  levels  1  and  3  in  the  model,  however,  is 
dominated  by  the  latent  heat  released  during  the  simulated  condensa¬ 
tion  (or  precipitation) .  This  process  accounts  for  some  70  percent 
of  the  heating  (as  opposed  to  cooling  in  the  model) ,  and  its  global 
distribution  may  be  obtained  from  the  precipitation  rates  of  Fig.  37 
by  conversion  to  the  equivalent  latent  heating  rate.  The  other  heat¬ 
ing  processes  important  in  the  model’s  tropospheric  heat  budget  are 
the  net  absorption  of  solar  insolation  and  the  net  flux  of  longwave 
radiation.  Lacking  appropriate  global  observations  of  these  processes, 
a  comparison  of  the  simulated  zonal  averages  with  the  available  ob¬ 
served  zonal  averages  has  been  given  elsewhere  (Gates,  1972).  A  more 
comprehensive  analysis  of  the  model’s  simulated  radiation  budget  is  in 
preparation  (Kahle,  1972). 

SURFACE  SENSIBLE-HEAT  FLUX 

The  upward  (turbulent)  flux  of  sensible  heat  at  the  surface  is 
given  by  the  expression 


r  =  p.Cnc  W(T  -  T.)  (14) 

K4  D  p  g  4 

where  T^  is  the  ground  surface  temperature,  T^  is  the  air  temperature 
just  above  the  ground,  and  p^,  C^,  and  W  are  as  in  (13).  Over  the 
oceans  T^  is  the  assigned  sea-surface  temperature  (see  Fig.  2),  whereas 
over  land  it  is  computed  from  an  assumed  condition  of  zero  net  heat 
flux.  The  surface  air  temperature  T^  resembles  the  model’s  free  air 
temperature  extrapolated  to  the  surface  but  is  also  affected  by  the 
simulated  low-level  convective  heat  transport.  The  sensible  heating 
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given  by  (14)  is  assumed  to  represent  a  flux  from  the  surface  to  the 
model’s  lower  level  (level  3),  where  it  serves  to  change  the  free- air 
temperature,  ,  as  a  part  of  the  diabatic  heating,  H^,  in  (3). 

The  simulated  surface  sensible-heat  flux  for  January  is  shown  in 
Fig.  44,  along  with  the  observed  pattern  as  given  by  Schutz  and  Gates 
(1971)  based  on  data  of  Budyko  (1963).  The  prominent  features  of  the 
simulation  (Fig.  44a)  are  the  maxima  of  upward  heat  flux  off  the  east 
coasts  of  North  America  and  Asia,  where  there  is  a  characteristic  flow 
of  relatively  cold  air  from  the  continents  over  the  warmer  oceans. 

While  these  features  are  realistic  in  both  pattern  and  position,  they 
are  approximately  twice  as  intense  as  the  fluxes  derived  from  the  ob¬ 
served  data  (Fig.  44b). 

A  second  noteworthy  feature  of  the  simulated  sensible  surface 
heating  is  the  calculation  of  a  negative  flux  over  much  of  the  tropical 
oceans,  with  values  as  large  as  200  ly  day  ^  off  the  central  west  coasts 
of  Africa  and  South  America.  According  to  Budyko’ s  (1963)  data  as  shown 
in  Fig.  44b,  the  observed  heat  flux  is  here  small  but  positive.  The 
overall  pattern  of  the  simulated  sensible— heat  flux  over  the  oceans 
(Fig.  44a)  bears  a  marked  resemblance  to  the  simulated  cloudiness 
(Fig.  34);  both  fields  are  dominated  by  low-level  convection  over  the 
sea.  Over  the  central  and  southern  parts  of  South  America  and  Africa, 
as  well  as  over  Australia,  the  model  has  successfully  simulated  upward 
heat  fluxes ,  but  here  again  the  maxima  are  approximately  twice  those 
observed  over  these  continents  in  the  southern  summer. 

The  errors  of  the  simulation  of  sensible— heat  flux  are  shown  more 
clearly  in  Fig.  45,  where  local  differences  of  more  than  100  ly  day 
are  not  uncommon.  In  Fig.  46,  the  simulated  flux  is  seen  to  average 
about  50  ly  day  ^  too  low  between  about  20°N  and  20°S,  and  is  more 
than  100  ly  day  ^  below  that  observed  in  high  northern  latitudes.  In 
contrast  to  the  observed  zonal  average  heat  flux  which  is  everywhere 
positive,  the  model  yields  a  downward  heat  flux  (from  the  air  to  the 
ground)  in  both  low  and  high  latitudes.  The  average  downward  (negative) 
fluxes  simulated  at  latitudes  north  of  about  50°N  are  probably  due  to 
the  model’s  calculation  of  low  surface  temperatures  over  the  continents, 
whereas  Budyko ’s  data  in  this  region  are  almost  entirely  over  the 
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oceans  (see  Fig,  44b),  If  we  compare  the  simulated  surface  heat  flux 
as  zonally  averaged  only  over  the  oceans  (shown  by  the  open  circles  in 
Fig.  46),  we  find  at  least  agreement  in  sign  with  Budykofs  data  at 
latitudes  north  of  50°N. 

It  appears  that  the  low  (and  even  negative)  zonally  averaged 
sensible-heat  flux  simulated  between  about  20°N  and  20°S  (Fig.  46)  is 
due  to  the  model's  systematic  overprediction  of  the  surface  air  tem¬ 
perature  over  the  tropical  oceans.  A  comparison  of  Figs.  3  and  13 
shows  that  the  simulated  surface  air  temperature  not  only  is  generally 
higher  than  that  observed  but  is  several  degrees  higher  than  the  pre¬ 
scribed  sea-surface  temperature  itself  over  all  tropical  oceans  between 
about  20°N  and  20°S.  This  has  resulted  through  (14)  in  the  reversal  of 
the  normal  (observed)  upward  sensible-heat  flux  in  this  region,  as 
shown  by  the  oceanic  zonal  averages  (open  circles)  in  Fig.  46.  The  use 
of  more  realistic  sea— surface  temperatures  and  a  more  detailed  treat¬ 
ment  of  the  surface  layer  in  the  model  could  be  expected  to  reduce 
these  discrepancies . 

SURFACE  SHORTWAVE  RADIATION 

Depending  on  wave  length,  the  shortwave  (solar)  radiation  incident 
at  the  top  of  the  model  atmosphere  is  either  subject  to  atmospheric 
absorption  or  subject  to  scattering,  both  by  the  air  and  by  clouds. 

Both  components  of  the  radiation  are  also  partially  reflected  at  the 
earth's  surface,  depending  upon  the  locally  prescribed  albedo,  with 
the  remainder  being  absorbed  at  the  surface.  The  global  distribution 
of  the  solar  radiation  incident  at  the  surface,  as  simulated  in  the 
model  for  January,  is  shown  in  Fig.  47a.  The  observed  distribution 
of  the  solar  radiation  received  at  the  surface  in  January  is  shown  in 
Fig.  47b,  based  on  the  data  of  Budyko  (1963)  as  summarized  by  Schutz 
and  Gates  (1971) . 

Both  the  observed  and  simulated  surface  solar  radiation  show  maxima 
over  the  southern  hemisphere  oceans,  as  expected  for  January,  with  a 
more-or-less  uniform  northward  decrease.  In  general,  however,  the 
simulated  surface  solar  radiation  (Fig.  47a)  is  somewhat  greater  than 
the  observed  incident  radiation  (Fig,  47b)  ,  as  may  be  more  clearly  seen 


Fig.  44a 


The  simulated  sensible  heat  flux  at  the  surface  for  January,  with  isolines  at  inter 
vals  of  100  ly  day--*-  and  the  zero  isoline  dashed. 


Fig.  44b  —  The  observed  sensible 
intervals  and  the  zer 
tional  data  are  missi 


January,  with  isolines  a 
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Fig.  46  The_ zonal  average  of  the  surface  sensible  heat  flux  as  simulated  for  January 
(full  line)  and  as  observed  (dashed  line)  according  to  the  data  of  Fig.  44. 
open  circles  denote  the  zonal  average  over  the  oceans  only. 


The 


-■88- 


in  the  difference,  or  simulation  error,  in  Fig.  48.  On  the  basis  of 
the  zonal  averages  shown  in  Fig.  49,  this  oversimulation  of  the  short¬ 
wave,  or  solar,  radiation  incident  at  the  surface  amounts  to  some 
134  ly  day  ^  in  the  global  mean.  Since  the  global  average  observed 
(between  «70°N  and  62®S)  is  377  ly  day  \  this  error  is  a  substantial 
fraction  of  the  solar  radiation  budget.  The  solar  radiation  absorbed 
at  the  surface  may  be  found  from  the  simulated  incident  radiation  data 
in  Figs.  47a  and  48  by  multiplying  by  (1  -  a) ,  where  a  is  the  surface 
albedo  (0.07  over  ocean,  0.14  over  land,  and  varying  with  latitude  from 
about  0.4  to  0.9  over  ice;  see,  however,  the  footnote  on  p.  8). 

This  evidence  suggests  that  the  model  permits  too  much  shortwave 
radiation  to  penetrate  to  the  surface  by  the  simulation  of  too  little 
absorption  or  reflection  in  the  troposphere.  Although  the  model’s 
simulated  radiation  budget  is  not  analyzed  in  detail  here,  the  analysis 
of  a  closely  related  experiment  has  shown  that  the  total  shortwave 
radiation  absorbed  in  the  atmospheric  column  is  reasonably  accurate 
(Gates,  1972).  Since  this  same  analysis  also  showed  that  the  simulated 
net  longwave  cooling  of  the  (tropospheric)  column  was  reasonably  close 
to  observation,  we  are  led  to  speculate  that  the  model  significantly 
underestimates  the  shortwave  radiation  reflected  to  space  by  clouds  . 

The  simulated  total  cloudiness  was  shown  earlier  (Figs.  34  to  36)  to 
be  systematically  too  low  by  about  0.2  coverage.  Noting  that  this 
error  is  half  the  observed  global  average  cloudiness  of  0.4,  the  model’s 
clouds  thus  reflect  on  the  average  only  about  half  as  much  solar  radia¬ 
tion  to  space  as  do  the  observed  clouds.  Since  reflection  by  clouds 
is  an  important  part  of  the  atmosphere’s  shortwave  radiation  budget, 
an  increase  in  the  model’s  cloudiness  by  about  a  factor  of  two  could 
therefore  be  expected  to  substantially  reduce  the  surface  shortwave 
radiation  error. 

SURFACE  HEAT  BALANCE 

Letting  the  solar  radiation  absorbed  at  the  surface  be  denoted  by 

S  ,  and  the  net  longwave  radiation  emitted  from  the  surface  by  R^,  the 
§ 

net  surface  heating  or  heat  balance  is  given  by  the  expression 

S  -  R,  -  T  -  Lp  E 
g  4  w 


(15) 
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where  T  is  the  surface  sensible-heat  flux  as  in  (14)  ,  E  is  the  evapor¬ 
ation  rate  as  in  (13),  and  L  is  the  latent  heat  of  condensation.  Over 
land  areas  in  the  model  it  is  assumed  that  (15)  is  equal  to  zero,  i.e., 
that  there  is  no  net  heating  at  the  ground  surface.  The  resulting 
equation  is  in  fact  used  to  determine  the  unknown  ground  temperature, 

Tg,  upon  which  R^  and  V  depend.  Over  the  oceans,  on  the  other  hand, 
no  such  balance  is  assumed  and  the  water's  surface  temperature  is  pre¬ 
scribed  (as  in  Fig.  3). 

The  simulated  surface  heating  (15)  for  January  is  shown  in  Fig.  50, 
along  with  the  observed  distribution  given  by  Budyko  (1963)  ,  as  sum¬ 
marized  by  Schutz  and  Gates  (1971) .  The  principal  areas  of  observed 
net  surface  heating  in  the  western  North  Atlantic  and  western  North 
Pacific  oceans  are  reasonably  well  simulated  by  the  model,  as  are  the 
regions  of  net  heat  loss  in  the  lower  latitudes.  Over  the  southern 
oceans  (say  south  of  40°S)  the  model  simulates  a  net  heating  several 
times  greater  than  that  given  by  the  data  of  Fig.  50b ,  as  is  readily 
seen  in  the  differences  between  the  simulated  and  observed  heat-balance 
distributions  given  in  Fig.  51.  A  local  error  of  300  ly  day  is  here 
not  uncommon,  and  is  in  many  instances  comparable  to  the  individual 
terms  of  (15) .  Over  the  continents  both  the  observed  and  simulated 
net  surface  heating  are  generally  slight.  If  we  recall  that  the  model 
was  designed  with  a  zero  surface  heat  balance  over  land,  the  small 
values  (<  10  ly  day  ^)  in  Fig.  50b  over  the  continents  are  thus  seen 
to  be  due  to  numerical  error. 

In  the  zonal  averages  shown  in  Fig.  52,  the  simulated  surface  heat 
balance  can  be  considered  a  good  approximation  to  observation  every¬ 
where  north  of  about  45°S.  In  diagnosing  this  balance  in  terms  of  (15), 
the  surface  longwave  radiation,  R^,  may  be  considered  to  be  modeled  with 
good  accuracy  at  all  except  polar  latitudes  in  view  of  the  generally 
successful  surface  temperature  simulation  (see  Figs.  13  to  15).  In 
the  tropics  the  excessive  solar  radiation  absorbed  at  the  surface,  S^, 
and  the  positive  heating  error  introduced  by  the  simulated  (negative) 
sensible-heat  flux,  F,  are  approximately  offset  by  the  excessive  heat, 
Lp^E,  lost  through  too  high  a  simulated  surface  evaporation  (see 
Fig.  42) .  In  the  middle  northern  latitudes  an  approximately  correct 


Fig.  47b  —  The  observed  shortwave  (solar)  radiation  incident  at  the  surface  in  January,  with 
isolines  every  25  ly  day” ^  and  the  250  ly  day  ^  isoline  dashed.  From  data  of 
Budyko  (1963)  as  summarized  by  Schutz  and  Gates  (1971) . 
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Fig.  49  —  The  zonal  average  of  the  simulated  shortwave  radiation  incident  at  the  surface 
(full  line)  and  that  observed  (dashed  line),  based  on  the  data  of  Fig.  47. 


Fig.  50a  —  The  simulated  surface  heat  balance  in  January,  given  by  the  excess  of  the  solar  radiation 
absorbed  over  the  (usual)  losses  due  to  evaporative  cooling,  long-wave  radiation,  and  up¬ 
ward  sensible  heat  flux*  The  isoline  interval  is  100  ly  day~l,  with  the  zero  isoline 
dashed.  Data  from  Figs.  40,  44,  and  47. 


The  observed  surface  heat  balance  in  January,  given  by  the  excess  of  the  solar  radiation 
absorbed  over  the  (usual)  losses  due  to  evaporative  cooling,  long-wave  radiation,  and  up 
ward  sensible  heat  flux.  The  isoline  interval  is  100  ly  day-1,  with  the  zero  isoline 
dashed.  Based  on  data  of  Budyko  (1963)  as  summarized  by  Schutz  and  Gates  (1971). 
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ace  heat  balance  as  simulated  for  January  (full 
line),  based  on  the  data  of  Fig.  50. 
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surface  heat  balance  results  from  compensation  between  a  simulated  T 
which  is  too  small  and  a  simulated  evaporation  which  is  too  great. 

Over  the  southern  hemisphere  oceans  (say  south  of  45°S)  ,  such 
error  compensation  does  not  occur.  The  excessive  positive  net  surface 
heating  simulated  between  about  50°S  and  70°S  (see  Figs.  50a  and  52) 
is  a  direct  result  of  the  simulation  of  excessively  great  solar-radiation 
absorption  at  the  surface,  as  may  be  seen  in  Fig.  49.  This  error  is 
evidently  due  to  the  small  amounts  of  simulated  cloudiness,  and  not  to 
the  inadvertent  surface-albedo  error  noted  earlier.  This  excessive 
solar  radiation  incident  (and  absorbed)  at  the  surface  is  responsible 
for  the  high  simulated  global  average  surface  heating  of  35  ly  day 
The  corresponding  observed  value  from  the  (incomplete)  global  data  of 
Fig.  50b  is  7  ly  day 
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VII.  SUMMARY  AND  CONCLUSIONS 


In  this  report,  the  January  global  climate  simulated  by  the  two- 
level  Mint  z-Arakawa  model  has  been  presented,  both  in  terms  of  the 
mean  global  patterns  and  the  associated  zonal  averages.  In  addition 
to  the  familiar  meteorological  variables  of  pressure,  temperature, 
wind,  humidity,  and  precipitation,  the  simulated  patterns  of  certain 
elements  of  the  moisture  and  energy  balances  are  also  given.  In  gen¬ 
eral,  those  fields  selected  were  ones  for.  which  global  observations 
were  available.  To  facilitate  comparison,  both  the  model’s  simulations 
and  the  observed  data  have  been  presented  in  a  uniform  format  (Schutz 
and  Gates,  1971);  the  difference  between  the  two  fields  (the  simulation 
error)  is  also  given. 

Considering  the  relative  simplicity  of  the  model,  the  simulated 

climate  may,  in  general,  be  considered  reasonably  accurate.  A  number 

of  characteristics  of  the  errors  of  the  January  simulation  are  sum- 

* 

marized  in  Table  1.  A  detailed  discussion  of  each  of  the  simulation 
errors  is  given  in  the  text,  along  with  the  complete  global  and  zonal 
average  error  distributions.  The  text  should  also  be  consulted  for  a 
discussion  of  the  simulated  mean  meridional  circulation  and  the  surface 
absorption  of  solar  radiation. 

Each  of  the  simulated  fields  shows  a  large-scale  pattern  of  de¬ 
parture  from  the  observed  distribution,  and  in  some  variables  this 
error  represents  a  substantial  fraction  of  the  observed  value  itself. 
This  is  particularly  true  for  zonal  wind,  relative  humidity,  precipita¬ 
tion,  evaporation,  and  sensible-heat  flux  and  heat  balance  at  the  sur¬ 
face.  Most  of  these  errors  are  systematic  with  respect  to  latitude 
and  may  be  seen  in  the  zonal  averages.  For  example,  the  simulated 
temperature  at  400  mb  is  about  6  deg  C  too  high  in  the  tropics  (Fig. 
21);  the  maximum  zonal  winds  (jet  stream)  in  the  northern  hemisphere 

& 

Excluded  from  this  summary  are  the  latitudes  south  of  70° S,  in 
view  of  the  inadvertent  ice-albedo  error  noted  earlier.  As  is  de¬ 
scribed  previously,  this  error  has  introduced  an  anomaly  of  10  to  15 
deg  C  in  the  surface  temperatures  over  Antarctica,  and  has  resulted 
in  the  simulation  of  excessive  evaporation  and  (convective)  precipita¬ 
tion  between  about  70°S  and  90° S. 
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Table  1 

SUMMARY  OF  ERRORS  OF  JANUARY  CLIMATE  SIMULATED  BY  THE  TWO-LEVEL  MINTZ-ARAKAWA  MODEL 


Error  Characteristics 
(Simulation  Minus  Observation) 

Variable 

Typical 

Local 

Error 

Magnitude 

Maximum 

Zonal  Average 
Error3  (and 
Simulated  Value) 

Global  Average 
Error  (and 
Simulated  Value) 

Remarks  on  Simulation 

Sea- level  pressure  (mb) 

Geopotential  height  (m) 

±10 

-10.3 

(1003.5) 

-0.8b 

(1010.3) 

midlatitude  quasi-statlonary  systems  too 
intense;  transient  cyclones  too  weak. 

800  mb 

+300  [ 

+382 

(2042) 

+300 

(2268)  \ 

absolute  heights  everywhere  too  high  ; 

400  mb 

+300 

+441 

(7795) 

+323 

(7600)  j 

midlatitude  gradient  too  steep. 

Temperature:  (deg  C) 

surface 

±6 

+4.0 

(-9.5) 

+1.2 

(13.4) 

tropical  oceans  and  continental  interiors 
too  warm. 

800  mb 

±6 

-5.3 

(2.3) 

-0.9 

(3.2) 

oceanic  regions  too  cold;  midlatitude 
gradient  too  great. 

400  mb 

±10 

+9.4 

(-18.5) 

+3.7 

(-24.2) 

tropics  too  warm;  midlatitude  gradient 
too  great. 

Zonal  wind  (m  sec 

800  mb 

400  mb 

Meridional  circulation 
(1012  g  sec-1) 

±10 

±20 

±40® 

+6.5 

+14.9 

-99 

(13.7) 

(35.1) 

(-101) 

+1*  °j 

+3.1d 

(3.1)  1 
(11.8)1 

f 

midlatitude  westerlies  too  strong. 

southern  hemisphere  Hadley  and  Ferrel 
cells  too  strong;  northern  Ferrel  cell 
too  weak. 

Relative  tumidity  (Z) 

tropics  too  dry,  midlatitudes  too  moist. 

800  mb 

±30 

-27.0 

(31.6) 

+1.4 

(59.2) 

Cloudiness 

-0.3 

-0.37 

(0.06) 

-0.26f 

(0.26) 

total  cloudiness  everywhere  too  little. 

Precipitation  (mm  day  1) 

±5 

+6. 38 

(7.4) 

+1.9® 

(4.1) 

tropical  rainfall  maxima  too  high  and 
too  far  poleward. 

Evaporation  (mm  day 

±5 

+3.0 

(6.8) 

+1.6h 

(4.0) 

excessive  evaporation  over  tropical 
oceans. 

Surface  sensible-heat 

flux  (ly  day~l) 

Surface  shortwave 

±100 

-124 

(-89) 

-17h 

(15) 

negative  flux  over  tropical  oceans;  up¬ 
ward  flux  over  continents  and  higher- 
latitude  oceans  too  great. 

radiation  (ly/day  ) 

±200 

+386 

(704) 

+1341 

(487) 

excessive  insolation  in  tropics  and  over 
southern  oceans. 

Surface  heat.,  balance 

(ly  day”1) 

±300 

+355 

(548) 

+17h 

(35) 

excessive  net  heating  over  southern 
oceans. 

Excluding  Antarctica. 

Due  to  arbitrary  reduction  to  sea  level. 
CSee  pp.  18  and  19. 
d Average  between  4*N  and  88°N. 
eAverage  between  60° S  and  60  *N. 


^Average  between  18°N  and  90°N. 
®See  footnote,  p.  67. 

^Average  between  58° S  and  78®N. 
* Average  between  62  °S  and  70°N. 
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are  about  60  percent  stronger  than  the  observed  winds  (with  a  smaller 
error  in  the  southern  (summer)  hemisphere,  Figs*  24  and  25);  the  rela¬ 
tive  humidity  is  about  half  that  observed  in  the  tropics  and  about  50 
percent  greater  than  the  observed  values  in  middle  and  high  latitudes 
(Fig,  33);  both  the  precipitation  and  the  evaporation  are  about  twice 
those  observed  in  the  tropics  (Figs,  39  and  42) ;  and  the  simulated  sur¬ 
face  sensible-heat  flux  is  negative  at  low  and  high  latitudes,  whereas 
the  observed  flux  is  everywhere  positive  (Fig.  46). 

Aside  from  the  systematic  bias  in  the  geopotential  noted  earlier, 
only  the  simulated  cloudiness  shows  an  error  of  uniform  sign;  judged 
from  the  best  available  observations  (U.S.  Air  Force;  see  Fig.  36), 
the  total  cloud  cover  simulated  by  the  model  is  only  about  half  that 
observed,  with  zonal  averages  below  the  observed  values  everywhere  in 
the  northern  hemisphere. 

From  such  an  overview  of  the  simulation  errors,  it  is  tempting  to 
speculate  on  their  causes  and  possible  interrelations.  On  the  basis 
of  the  present  evidence,  there  appear  to  be  two  major  sources  of  signif¬ 
icant  error  in  the  model.  First  are  the  low  values  of  simulated  cloud¬ 
iness  just  mentioned.  The  modelrs  cloudiness  is  primarily  determined 
by  the  simulated  occurrence  of  low-level  convection,  which  occurs  almost 
exclusively  over  the  oceans  (see  Fig.  36).  Although  no  separate  low 
cloudiness  observations  are  available,  it  is  possible  that  this  cloud 
type  is  reasonably  well  accounted  for  in  the  model.  The  substantial 
amounts  of  middle-  and  high-level  clouds  observed  in  the  atmosphere, 
however,  are  represented  in  the  model  by  relatively  small  cloud  amounts, 
which,  moreover,  are  simulated  only  during  condensation.  The  present 
model  makes  no  attempt  to  simulate  nonprecipitating  clouds,  and  it  is 
the  neglect  of  this  cloudiness  which  is  likely  responsible  for  the  low 
total  cloud  cover  simulated.  According  to  London  (1957)  the  energy 
reflected  from  clouds  is  approximately  half  the  shortwave  radiation 
absorbed  at  the  surface.  The  model's  low  simulated  cloudiness  thus 
reflects  too  little  of  the  incoming  solar  radiation  and,  thereby,  per¬ 
mits  too  much  radiation  to  be  absorbed  at  the  surface  (Fig.  49).  In 
January  this  excessive  solar  radiation  absorption  is  a  maximum  in  the 
southern  hemisphere  and  is  responsibe  for  the  most  prominent  error  of 
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the  simulated  surface  heat  budget  (Fig.  52).  Most  of  this  energy, 
however,  disappears  into  the  ocean,  which  is  an  infinite  heat  reservoir 
in  the  model. 

A  second  significant  error  in  the  model  is  the  high  values  of  pre¬ 
cipitation  simulated  in  the  tropical  region  between  about  25 °N  and  25° S. 
As  was  previously  noted,  almost  all  of  this  precipitation  is  determined 
by  the  simulated  occurrence  of  penetrating  convection.  This  excessive 
rainfall  is  accompanied  by  the  release  of  too  much  latent  heat  of  con¬ 
densation,  and  may  therefore  be  responsible  for  the  excessively  warm 
air  simulated  at  400  mb  (Fig.  21).  The  midlatitude  meridional  temper¬ 
ature  gradient  is  thereby  increased,  which  in  turn  contributes  to  the 
simulation  of  excessively  strong  westerly  winds  and  too  steep  a  merid¬ 
ional  geopotential  height  gradient.  Rather  than  increasing  the  insta¬ 
bility  of  the  transient  baroclinic  waves  in  middle  latitudes  as  might 
be  expected,  these  errors  appear  to  have  intensified  the  quasi-s tat ionary 
long  waves  and  resulted  in  an  excessively  strong  mean  meridional  circu¬ 
lation  (Fig.  30b).  Although  it  is  dangerous  to  speculate  in  a  qualita¬ 
tive  manner,  it  may  be  that  the  excessive  tropical  precipitation  removes 
too  much  moisture  from  the  atmosphere  (modeled  only  at  the  lower  level) 
and  is  thus  a  cause  of  the  low  simulated  relative  humidity  at  800  mb 
between  about  25°N  and  25°S  (Fig.  33).  This  in  turn  would  promote  the 
excessive  evaporation  in  the  same  region  (Fig.  42),  since  the  surface 
moisture  is  fixed  over  the  oceans  where  nearly  all  of  the  evaporation 
occurs. 

If  the  major  errors  of  the  present  control  integration  are  indeed 
associated  with  the  model’s  underprediction  of  the  cloudiness  and  its 
overprediction  of  the  convection  as  hypothesized  above,  then  the  appro¬ 
priate  revision  of  these  parameterizations  could  be  expected  to  markedly 
improve  the  accuracy  of  the  climatic  simulation.  Accordingly,  research 
is  now  under  way  to  increase  the  amount  of  cloudiness  simulated  in  the 
model  by  suitable  revision  of  the  relation  (12)  and  by  the  introduction 
of  a  provision  to  allow  midlevel  nonconvective  cloudiness  even  in  the 
absence  of  precipitation.  Other  research  is  in  progress  whose  aim  is 
to  reduce  the  amount  of  convective  precipitation  by  suitable  revision 
of  the  convective  parameterization.  (This  parameterization  has  also 
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been  revised  in  the  design  of  a  new  three-level  version  of  this  model 

[Arakawa,  1972].)  With  these  modifications,  a  new  January  integration 

* 

is  envisaged,  which  will  also  include  the  appropriate  monthly  distribu¬ 
tion  of  sea-surface  temperature  (replacing  the  annual  averages  used 
here) ,  and  a  more  realistic  distribution  of  the  seasonal  surface  albedo 
(replacing  the  values  0.14  and  0.07  now  assigned  to  all  land  and  ocean, 
respectively,  and  the  (erroneous)  values  used  for  snow  and  ice) . 

In  interpreting  the  results  of  this  or  any  other  general  circulation 
model,  it  is  important  to  recognize  that  there  may  be  significant  varia¬ 
tions  in  the  simulated  seasonal  averages  from  one  year  to  the  next;  most 
analyses,  including  the  present  one,  have  considered  only  a  single  month’s 
simulation.  Some  information  on  annual  variability  of  the  present  model 
over  a  three-year  period  is  in  preparation  by  Mintz  et  al.  (1972),  and 
longer  runs,  are  envisaged  with  improved  versions  of  the  model  as  dis¬ 
cussed  above.  Further  work  of  this  sort  is  needed  in  order  to  determine 
more  accurately  the  model’s  true  ability  to  simulate  climate  and,  there¬ 
by,  to  assess  its  usefulness  as  a  tool  in  the  study  of  climatic  change. 
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